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Abstract. This paper concerns the hnk between the dynamical behaviour of a dynam- 
ical system and the dynamical behaviour of its numerical simulations. Here, we model 
numerical truncation as a spatial discretization of the system. Some previous works 
on well chosen examples (see e.g. [Boy86], [GT83], [GB88], [KMP99] or [Lan98]) show 
that the dynamical behaviours of dynamical systems and of their discretizations can be 
quite different. In this paper we are interested in generic homeomorphisms of compact 
manifolds. So our aim is to tackle the following question: can the dynamical properties 
of a generic homeomorphism be detected on the spatial discretizations of this homeo- 
morphism? We study both conservative and dissipative homeomorphisms (i.e. with or 
without the assumption that the homeomorphism preserves a given measure). 

We will prove that the dynamics of a single discretization of a generic conservative 
homeomorphism does not depend on the homeomorphism itself, but rather on the grid 
used for the discretization. Therefore, dynamical properties of a given generic conserva- 
tive homeomorphism cannot be detected using a single discretization. Nevertheless, we 
will also prove that some dynamical features of a generic conservative homeomorphism 
(such as the set of the periods of all periodic points) can be read on a sequence of 
finer and finer discretizations. The dissipative case is somehow opposite: we will see 
that the dynamics of a generic dissipative homeomorphism is well-approximated by the 
dynamics of any fine enough discretization of this homeomorphism. 
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1. Introduction 

When one tries to simulate a discrete time dynamical system with a computer, calcu- 
lations are performed with a finite numerical precision. Numerical errors made at each 
iteration may add up, so that after a while, the numerically calculated orbit of a point will 
have nothing in common with the actual one. Nevertheless, a numerically calculated orbit 
is close to a actual orbit at any time, thus one can hope that the collective behaviour of 
numerically calculated orbits provides information about the collective behaviour of actual 
orbits. Interesting properties of dynamical systems are often asymptotic (i.e. concern the 
behaviour of the system in the long term) on the one hand, and collective (i.e. concern the 
collective behaviour of orbits or even the action of the system on open sets and probability 
measures) on the other. The purpose of our study is to tackle the question: can such 
properties be detected on numerical simulations? 

We will model numerical truncation as a spatial discretization^. Consider a dynamical 
system whose phase space is a manifold X and whose evolution law is given by a map / 
from X into itself. Roughly speaking, a numerical simulation of the system with a precision 
IQ-N replaces the continuous phase space X by a discrete space En made of points of X 
whose coordinates are decimal numbers with at most N decimal places, and replaces the 
map / by its discretization /n : En ^ En that maps x e En to the point (or one of 
the points) of En nearest f{x). We can clarify a bit the question raised in the previous 
paragraph: 

Question. Can the dynamical properties of a homeomorphism / be read on the dynamics 
of its discretizations (/Ar)Ar>o? 

Questions about discretizations of dynamical systems have mainly been studied em- 
pirically, using numerical simulations of well-chosen examples (see e.g. [Boy86], [GB88], 
[KMP99] or [Lan98]). In particular, a short article of J.-M. Gambaudo and C. Tresser 
[GT83] gives the examples of two simple homeomorphisms of [0, 1]^, which both have sinks 
that attract most of the orbits. However, those sinks are undetectable in practice, simply 
because the connected components of their basins are much too small. Thus, the actual 
dynamics of these homeomorphisms does not have much in common with that observed 
on simulations. From a more theoretical point of view, E. Ghys studied (in [Ghy94]) the 
dynamical behaviour of discretizations of the Anosov automorphism of the torus given 
by (x^y) {y^x -\-y). E. Ghys notes that the dynamics of these discretizations (he shows 
that they are permutations with remarkably small orders) has nothing to do with the actual 
dynamics of the automorphism, which is a paradigm of chaotic dynamics. The very specific 
arithmetical properties of this example let one hope that such a difference between the 
actual dynamics of a system and the dynamics of its discretizations might be exceptional. 
This hope is supported by more recent results of S. Hayashi in [Hayl2], roughly speaking, 
he shows that every C^-diffeomorphism of a compact manifold with no positive Lyapunov 
exponent can be approximated by another diffeomorphism whose sinks are "observable in 
practise". Furthermore, T. Miernowski made a fairly comprehensive study of discretizations 
of circle homeomorphisms [Mie06]. He basically shows that the dynamical properties of a 
typical (generic or prevalent) circle homeomorphism/diffeomorphism can be read on the 
dynamics of its discretizations. All these results suggest to tackle the above question not 
for arbitrary homeomorphisms but rather for typical homeomorphisms. 

In this paper we study dynamical properties of discretizations of generic homeomor- 
phisms in the sense of Baire. We establish properties for both dissipative, i.e. arbitrary 
homeomorphisms of X, and conservative homeomorphisms, i.e. homeomorphisms of X 
that preserve a given good probability measure. 

The results we obtained concern generic homeomorphisms of a compact manifold (with 
boundary) of dimension n > 2. In this introduction we state the main results for the torus 
and for "uniform" discretizations on T"^. The general framework in which the results 
are valid is defined in section 2. 



See page 9 about the accuracy of this model. 
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Let n > 2. For all A/' G N, let En be the finite subset of the torus = R^/Z^ made 
of the points whose coordinates are decimal numbers with at most TV decimal places: 



Let Pn be a projection of T"^ on En minimizing the euclidian distance: for x G T^, Pn{x) 
is one of the points of En nearest x. Let /n - En En be the discretization of the 
homeomorphism / : according to En defined by /at = Pn o /• 

We will prove many results, concerning various aspects of the dynamics of the discretiza- 
tions. In the sequel of this introduction, we try to organize those results according to some 
"lessons": 

1) The dynamics of a single discretization of a generic homeomorphism has in general nothing 
to do with the dynamics of the initial homeomorphism. 

Given a finite set E^ the dynamics of any finite map a : E ^ E is quite simple: given 
X e E, the orbit {a^{x))k is preperiodic. Hencefore, the union of periodic orbits of cr, which 
we will denote by rt{a) and called the maximal invariant set of a, is exactely the union of 
the cj-limits sets of points of E. To study the dynamics of a one can focus on quantities 
such as the cardinal of ^^(cr), the stabilization time of a (i.e. the smallest t G N such that 
(J*{E) = r^(cr)), the number of orbits of (J\n(^a)^ their lengths, the period of (J\n(^a) - • • 

Another interesting dynamical quantity for a finite map a : E ^ E is the recurrence 
rate: it is the ratio between the cardinals of rt{a) and that of E. The starting point of our 
article is a question from E. Ghys (see chapter 6 of [Ghy94]): for a generic conservative 
homeomorphism of the torus, what is the asymptotical behaviour of the sequence of recur- 
rence rates of /at? A partial answer to this question was obtained by T. Miernowski in the 
second chapter of his thesis. 

Theorem (Miernowski). For a generic conservative homeomorphism'^ f , there are infin- 
itely many integers N such that the discretization /n is a cyclic permutation. 

To prove this theorem, T. Miernowski combines a genericity argument with a quite 
classical technique in generic dynamics of homeo(and auto)morphisms: approximation by 
permutations (see e.g. [KS70], [Alp 78], [Hal56]), and more precisely Lax's theorem (see 
[Lax71]), which states that any conservative homeomorphism of the torus can be approx- 
imated by cyclic permutations of the discretization grids. In fact this proof can be gen- 
eralized to obtain many results about discretizations. We will establish some variants of 
Lax theorem; each of them, combined with a genericity argument, leads to a result for 
discretizations of generic homeomorphisms. For instance: 

Theorem. For a generic conservative homeomorphism f , it exists C > such that there 
are infinitely many integers N such that the cardinal ofVt{fN) is smaller than C. 

Note that the combination of these two theorems answer the question of E. Ghys: for a 
generic homeomorphism /, the sequence of the recurrence rate of /at accumulates on both 
and 1, one can even show that it accumulates on the whole segment [0, 1]. 

Let us move on to the behaviour on the maximal invariant set (^(/at). Another variations 
of Lax's theorem lead to a theorem that enlightens the behaviour of the discretizations on 
their maximal invariant set: 

Theorem. For a generic conservative homeomorphism f , there are infinitely many integers 
N such that the discretization fN is a cyclic permutation and for all M G there are 
infinitely many integers N such that fN is a permutation of En having at least M periodic 
orbits. 



I.e. there is a Gs dense subset of the set of conservative homeomorphisms of the torus on which the 
conclusion of theorem holds. 
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To summarize, generically, infinitely many discretizations are cyclic permutations, but 
also infinitely many discretizations are highly non injective or else permutations with many 
cycles. In particular, it implies that for all x G X, it exists infinitely many integers N such 
the orbit of xn under fj^ does not shadow the orbit of x under /: in this sense, generically, 
the dynamics of discretizations does not reflect that of the homeomorphism. 

2) A dynamical property of a generic homeomorphism can not be deduced from the frequency 
it appears on discretizations either. 

The previous theorems express that the dynamics of a single discretization does not 
reflect the actual dynamics of the homeomorphism. However, one might reasonably expect 
that the properties of the homeomorphism are transmitted to many discretizations. More 
precisely, one may hope that given a property (P) about discretizations, if there are many 
N such that the discretization /n satisfies (P), then the homeomorphism satisfies a similary 
property. It is not so, for instance: 

Theorem. For a generic conservative homeomorphism f , the proportion of integers N 
between 1 and M such that /n is a cyclic permutation accumulates on both and 1 when 
M goes to infinity. 

In fact, for most of the properties considered in the previous paragraph, the frequency 
they appear on discretizations of orders smaller than M accumulates on both and 1 when 
M goes to infinity. 

3) Some important dynamical features of a generic homeomorphism can be detected by looking 
at some dynamical features of all the discretizations. 

We have observed that one can not detect the dynamics of a generic homeomorphism 
when looking at the dynamics of its discretizations, or even at the frequency some dynamics 
appears on discretizations. Nevertheless, some dynamical features can be deduced from 
the analogous dynamical features of all the discretizations. For instance, the periods of 
periodic orbits of a homeomorphism can be read on the periods of periodic orbits of its 
discretizations: 

Proposition. A generic homeomorphism f has a periodic orbit with period p if and only 
if there exists infinitely many integers N such that f^ has a periodic orbit with period p. 

We will also prove a theorem in the same vein for invariant measures. It expresses 
that the set of invariant measures of the homeomorphism can be deduced from the sets of 
invariant measures of its discretizations. More precisely: 

Theorem. Let Mn be the set of probability measures on En that are invariant under /at- 
For a generic conservative homeomorphism f , the upper limit over N (for Hausdorff topol- 
ogy) of the sets A4n is exactly the set of probability measures that are invariant under f . 

4) Physical measures of a generic homeomorphism can not be detected on discretizations. 

Given x G T^, the Birkhoff limit of x is defined (when it exists) as the limit of 
m ^T=Q fl^x when m goes to infinity; the basin of a measure /i is the set of points x 
whose Birkhoff limit coincides with /i. Heuristically, the basin of a measure is the set of 
points that can see the measure. A Borel measure is said physical if its basin has pos- 
itive Lebesgue measure (see e.g. [You02]). The heuristic idea underlying this concept is 
that physical measures are invariant measure which can be detected "experimentally" (since 
many initial conditions lead to these measures). 

According to this heuristic, one could expect physical measures to be the only invariant 
measures that can be detected on discretization. This is not the case: for a generic con- 
servative homeomorphism, there exists a unique physical measure, say Lebesgue measure. 
According to the previous theorem, invariant measures of the discretizations accumulate 
on all the invariant measures of the homeomorphism and not only on Lebesgue measure. 
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However, one can still hope to distinguish the physical measure from other invariant 
measures. For this purpose, we define the canonical physical measure associated to a 
discretization /at: it is the limit in the sense of Cesaro of the images of the uniform measure 
on En by the iterates of /at: if A at is the uniform measure on En, then 



This measure is supported by the maximal invariant set Q^/n)'', it is uniform on every 
periodic orbit and the total weight of a periodic orbit is proportional to the size of its 
basin of attraction. The following theorem expresses that these measures accumulate on 
the whole set of /-invariant measures: physical measures can not be distinguished from 
other invariant measures on discretizations, at least for generic homeomorphisms. 

Theorem. For a generic conservative homeomorphism f , the set of accumulations points 
of the sequence {/J^N)Ne'N is the set of all f -invariant measures. 

5) The dynamics of discretizations of a generic dissipative homeomorphism tends to that of the 
initial homeomorphism. 

We then study properties of discretizations of generic dissipative homeomorphisms^. 
The basic tool is the shredding lemma of F. Abdenur and M. Andersson [AA12], which 
implies that a generic homeomorphism has a "sink dynamics". This easily transmits to 
discretizations, for example the basins of attraction of the homeomorphism can be seen 
on all the fine enough discretizations. Moreover there is convergence of the dynamics of 
discretizations /at to that of /: 

Theorem. For a generic dissipative homeomorphism f , for all e > and all S > 0, there 
is a full measure dense open subset O of T^ such that for all x e O, all S > and all N 
large enough, the orbit of xn under /at S -shadows the orbit"^ of x under f . 

Thus, one can detect on discretizations the dynamics of a generic dissipative homeo- 
morphism, which is mainly characterized by position of the sinks and of the corresponding 
basins of attraction. Note that this behaviour is in the opposite of the conservative case, 
where the individual behaviour of discretizations does not indicate anything about the 
actual dynamics of the homeomorphism. 

6) In practice, one can not deduce the dynamics of a conservative homeomorphism from its 
discretizations, and the dynamics of a dissipative homeomorphism can be detected on discretiza- 
tions only if the basins are large enough. 

Finally, we compare our theoretical results with the reality of numerical simulations. 
Indeed, it is not clear that the behaviour predicted by our results can be observed on 
computable discretizations of an homeomorphism defined by a simple formula. On the one 
hand, all our results are valid "for generic homeomorphisms": nothing indicates that these 
results apply to practical examples of homeomorphisms defined by simple formulas. On the 
other hand, results such as "there are infinitely many integers N such that the discretization 
of order A/". . . " provide no control over the integers N involved, they may be so large that 
the associated discretizations are not calculable in practice. 

For conservative homeomorphisms, our numerical simulations produce mixed results. 
From a quantitative viewpoint, the behaviour predicted by our theoretical result cannot be 
observed on our numerical simulations. For example, we can not observe any discretization 
whose recurrence rate is equal to 1 (i.e. which is a permutation). This is nothing but 
surprising: the events pointed out by the theorems are a priori very rare. For instance, 
there is a very little proportion of bijective maps among maps from a given finite set into 




M-l 



m=0 



I.e. without assumption of preservation of a given measure. 
"^I.e. for all A; G N, d(f^(xN)J^ (x)) < 6. 
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itself. From a more qualitative viewpoint, the behaviour of the simulations is quite in ac- 
cordance with the predictions of the theorical results. For example, for a given conservative 
homeomorphism, the recurrence rate of a discretization does depends a lot on the size of 
the grid used for the discretization. Similarly, the canonical invariant measure associated 
with a discretization of a homeomorphism / does depend a lot on the size of the grid used 
for the discretization. 

We also carried out simulations of dissipative homeomorphisms. The results of dis- 
cretizations of a small perturbation of identity (in topology) may seem disappointing 
at first sight: the sinks of the initial homeomorphisms can not be detected, and there is 
little difference with the conservative case. This behaviour is similar to that highlighted 
by J.-M. Gambaudo and C. Tresser in [GT83] (see page 1). That is why it seemed to us 
useful to test an homeomorphism which is close to the identity, but whose basins are 
large enough. In this case the simulations point out a behaviour that is very similar to that 
described by theoretical results, namely that the dynamics converges to the dynamics of 
the initial homeomorphism. In fact, one have actually observed behaviours as described by 
theorems only for examples of homeomorphisms with a very few number of sinks. 

We carried out many other numerical simulations, they can be found on the web page: 
http : //www. math. u-psud. fr/~guiheneu/Simulat ions .html 

In section 2 we will present the framework. Sections 4 to 9 concern conservative homeo- 
morphisms. More precisely, the results concerning the behaviour of a single discretization 
are set out in section 6, and those concerning the average behaviour of discretizations in 
section 7. We then come to the results about all the discretizations in section 8 and about 
physical measures in section 9. The behaviour of discretizations of generic dissipative 
homeomorphisms is established in section 11. Finally, the results of numerical simulations 
are presented in part 3. 

2. Framework 

The manifold X and the measure A. The results stated in the introduction for the 
torus and the Lebesgue measure Leb actually extend to any smooth manifold X with 
dimension n > 2, compact and possibly with boundary, endowed with a Riemannian met- 
ric d. We fix once and for all such a manifold X endowed with the metric d. In the general 
case, Lebesgue measure on T"^ can be replaced by a good measure A on X: 

Definition 1. A Borel probability measure A on X is called a good measure^ or an Oxtoby- 
Ulam measure if it is nonatomic, it has total support (it is positive on nonempty open sets) 
and it is zero on the boundary of X. 

We fix once and for all a good measure X on X . 

Notation 2. We denote by Homeo(X) the set of homeomorphisms of X, endowed by the 
metric d defined by: 

d{f,g) = sup d{f{x),g{x)). 

xex 

We note Homeo(X, A) the subset of Homeo(X) made of the homeomorphisms that preserve 
the measure A, endowed with the same metric d. Elements of Homeo(X) will be called 
dissipative homeomorphisms and which of Homeo(X, A) conservative homeomorphisms. 

Generic properties in Homeo(X) and Homeo(X, A). The topological spaces Homeo(X) 
and Homeo(X, A) are Baire spaces (see [Guil2]), i.e.: the intersection of every countable 
collection of dense open sets is dense. We call Gs a countable intersection of open sets; a 
property satisfied on at least a Gs dense set is called generic. Note that in a Baire space, 
generic properties are stable under intersection. 

Sometimes we will use the phrase "for a generic homeomorphism / G Homeo(X) (resp. 
Homeo(X, A)), we have the property (P)". By that we will mean that "the property (P) is 
generic in Ilomeo(X) (resp. Homeo(X, A))", i.e. "it exists Gs dense subset G of Homeo(X) 
(resp. IIomeo(X, A)), such that every f ^ G satisfy the property (P)". 



DYNAMICAL PROPERTIES OF SPATIAL DISCRETIZATIONS OF A GENERIC HOMEOMORPHISM 7 



Discretization grids, discretizations of a homeomorphism. We now define a more 
general notion of discretization grid. Some of the assumptions about these grids will be 
useful later, it will be the subject of the next paragraph. 

Definition 3 (Discretization grids). A sequence of discretization grids on X is a sequence 
{EN)Ne'N of discrete subsets of X\dX^ such that the grids are more and more precise: for 
all £ > 0, it exists TVq G N such that for all N > Nq^ the grid En is e-dense. We denote by 
qN the cardinal of En- 

We fix once and for all a sequence {EN)Ne^ of discretization grids on X. We can now 
define discretizations associated to these grids: 

Notation 4 (Discretizations). Let Pat be a projection of X on En (the projection of 
xo G X on En is some G En minimizing the distance d{xo^y) when y runs through 
En)- Such a projection is uniquely defined out of the set made of the points x G X 
for which there exists at least two points minimizing the distance between x and En- On 
E^j^ the map Pn is chosen arbitrarily (nevertheless mesurably). For x e X we note xn 
the discretization of order N of x, defined hy xn = Pn{x). For / G Homeo(X) we note 
/at : En En the discretization of order X of /, defined by /at = Pn ^ /• 

Let Vn he the set of homeomorphisms g such that g{EN) H E'js^ = 0. Subsequently we 
will restrict the study to the Gs dense set Hatgn which is the set of homeomorphisms 
whose N-th discretization is uniquely defined for all X G N. 

If a : En En and / G Homeo(X), we note dN^f^cr) the distance between f\E^ and 
cr, considered as maps from En into X. 

Remark 5. One might ask why the points of the discretization grids are supposed to be 
inside X. The reason is simple: a homeomorphism f of X sends dX on dX. Putting points 
of some grids on the edge could perturb the dynamics of discretizations /at- In particular 
it would introduce at least one orbit^ with length smaller than Card(£^Ar H dX). 

Probability measures on X. From section 8, we will be interested in ergodic properties 
of discretizations of /. Note V the set of Borel probability measures on X endowed with 
the weak-star topology: a sequence (i^m)mGN of V tends to u G V (denoted by Um ^) if 
for all continuous function (/:) : X ^ R, 



Under these conditions the space V is metrizable and compact, therefore separable (Pro- 
horov, Banach-Alaoglu-Bourbaki theorem). 

To study ergodic properties of homeomorphisms and their discretizations, we define 
natural invariant probability measures associated with these maps: 

Definition 6. For any nonempty open subset U of X, we note Xu the normalized restriction 
of A on [/, i.e. Xu = j^j^'^lU' We also note Xn,u the uniform probability measure on En^U 
and Aat = AAr,^:- For x G X we note (when the limit exists) the Birkhoff limit of x: 





m—l 



and similarly for /at- When it is well defined, we set 




m—l 



and 




m—l 



Finally, we note ji^^ = ji^^ . 



'Recall that orbit means forward orbit. 
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We just define two types of invariant measures: on the one hand from a point x, the 
other from the uniform measure A. The link between it is done by the following proposition: 

Proposition 7. When U is an open set whose almost every point admit a Birkhoff limit, 
the measure /Jj(j is well defined and satisfies, for every continuous map cp : X ^ H, 



Similarly, 




Proof of proposition 7. It follows easily from the dominated convergence theorem. □ 

Hypothesis on discretization grids. Previously, we have given a very general definition 
of the concept of sequence of discretization grids. In some cases, we will need additional 
technical assumptions about these sequences of grids. Of course all of them will be satisfied 
by the uniform discretization grids on the torus (as defined in the introduction). 

The first assumption is useful to prove Lax's theorem (theorem 19), and therefore nec- 
essary only in the part concerning conservative homeomorphisms. 

Definition 8 (Well distributed and well ordered grids). We say that a sequence of dis- 
cretization grids (^Ar)ArGN is well distributed if one can associate to each x G En a subset 
Cn,x of ^1 which will be called a cube of order such that: 

• for all N and all x G En, x G Cn,x, 

• for all {CN,x}xeEN is a measurable partition of X: \Jj.^En ^^^x is full measure 
and for x, y two distinct points of En, the intersection Cn^x^Cn^v is null measure, 

• the diameter of the cubes of order N tends to 0: max^^^^;^ didim{CN,x) — > 0- 

If {EN)Ne'N is well distributed and if {CN,x}Ne'N,xeEN is a family of cubes as below, 
we will say that {EN)Ne'N is well ordered if, for a fixed the cubes {CN,x}xeEN can be 
indexed by Z/^atZ such that two consecutive cubes (in Z/^atZ) are adjacent. 

The two following definitions describe assumptions that will be useful especially to obtain 
properties in average. 

Definition 9 (Refining grids). We say that a sequence of discretization grids refines if for 
ah A^, TV' G N such that TV < TV', one has En C En^^ 

Definition 10 (Self-similar grids). We say that a sequence of discretization grids {EN)Ne'N 
is self similar if for all £ > 0, there exists TVq, TVi G N such that for all TV > TVi, the set En 
contains disjoint subsets Ej^, . . . , E^^ whose union fills a proportion greater than 1 — e of 
En, and such that for all j, Ej^ is the image of the grid Enq by a bijection hj which is 
5-close to identity. 

We say that a sequence of discretization grids {EN)Ne'N is strongly self similar if it is 
self similar and for all TV > TVq, one of the hj equals to identity. 

Remark 11. One easily verify that the hypothesis "being strongly self-similar" implies both 
"being self-similar" and "refining". 

3. Some examples of discretization grids 

In the previous section we set properties on discretizations — namely being well dis- 
tributed, being well ordered, refining, being strongly self-similar or self-similar — that will 
be used subsequently. In this section we give some examples of grids that verify some of 
these assumptions. 
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Figure 1. Uniform discretization grids of order 5 on the torus {E^ left) and on the 
cube P {E^ middle and Ej^ right) and their associated cubes 



Uniform discretization grids on the torus. The simplest example, which will be used 
for the simulations, is that of the torus T"^ = 11^ /Z^ of dimension n > 2 endowed with 
discretizations called uniform discretizations, defined from the fundamental domain = 
[0, 1]^ of T^: take an increasing sequence of integers (/cAr)ArGN and set 



E 



N 



e r\ Vj, <ij <kN - I 



c 



N,{ii/N,...,ir^/N) 



n 



kN 



2kN ' k 



N 



2kN 



We easily verify that this sequence of grids is well distributed, well ordered and self- 
similar. If we further assume that for any G N, /cat divides /cat+i (which is true when 
kN = with p > 2), then the sequence is strongly self-similar (therefore refines). When 
kN = with p = 2 (resp. p = 10) the discretization performs what one can expect from a 
numerical simulation: doing a binary (resp. decimal) discretization at order N is the same 
as truncating each binary (resp. decimal) coordinate of the point x e to the 7V-th digit, 
i.e. working with a fixed digital precision^. 

Uniform discretization grids on the cube. Similarly, one can define uniform discretiza- 
tions on the cube = [0, 1]"^ by setting 



< ij < kN 



1 



still with an increasing sequence {kN)Ne'N of integers, to which are associated the cubes'" 
(see figure 1) 



kN — 1 k 



N 



1 



^N,{h/N,...,ir,/N) = Y[ 

i=i 

As before, one easily verifies that this sequence of grids is well distributed, well ordered 
and self-similar (see figure 2). If we further assume that for any A/" G N, /cat divides /cat+i, 
then the sequence is strongly self-similar (therefore refines). 

One can also take discretizations according to the centers of the cubes. 

1/2 in + 1/2^ 



^N 



kN kN 
This time the cubes are much more natural (see figure 1): 



G Vj, 0<ij < k 



N 



c 



N,{ii/N,...,ir^/N) 



n 



kN 



^N 



Again, we easily verify that this sequence of grids is well distributed, well ordered and 
self-similar (but does not refine). 



In fact, in practice the computer works in floating point format, so that the number of decimal places 
is not the same when the number is close to or not. 

^Be careful, these cubes have their vertices on the grid of order k]\j — 1. 
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Figure 2. Self similarity of grids E% 

Discretization grids on an arbitrary manifold. In fact, discretizations Ep^ on can 
be generalized to an arbitrary manifold X with the Oxtoby-Ulam theorem (see [0U41]): 

Theorem 12 (Oxtoby, Ulam). Under the assumptions that have been made on X and X, 
it exists a map (j) : ^ X such that: 

(1) (j) is surjective 

(2) (j)^fn is an homeomorphism on its image 

(3) (j){dl^) is a closed subset of X with empty interior and disjoint from (/>(/"') 

(4) xmn) = 

(5) ^*(A) = Leb^ where Leb is Lebesgue measure. 

This theorem allows us to define discretization grids on X from uniform discretizations 
on the cube. We obtain the following informal proposition: 

Proposition 13. // there is a sequence {En)n of grids on whose elements are not on 
the edge of the cube, then its image by (j) defines a sequence of grids on X which satisfy the 
same properties as the initial grid on the cube. 

Remark 14. The example also includes the case where X = A = Leb and where the 
grids are the images of the grids E^ by a unique homeomorphism of X preserving Lebesgue 
measure. 

Part 1. Discretizations of a generic conservative homeomorphism 

Recall that we have fixed once and for all a manifold X, a good measure A on X and a 
sequence of discretization grids {EN)Ne'N (see definition 3). It is further assumed that this 
sequence of grids is well distributed and well ordered (see definition 8). In this section we 
will focus on discretizations of elements of Homeo(X, A), so homeomorphisms will always 
be supposed conservative. 

4. Dense types of approximations 

To begin with we define the notion of type of approximation. 

Definition 15. A type of approximation T = (T/v)ArGN is a sequence of subsets of the set 
T{En, En) made of applications from En into itself. 

Let U be an open subset of Homeo(X, A). A type of approximation T is said to be dense 
in Z// if for all / G Z//, 6 > and all Xq G N, there exists N > Nq and ctat G Tn such that 
dNif^cTN) < s (recall that c^at is the distance between f\EN ^^a'^ considered as maps 
from En into X). 

The goal of this paragraph is to obtain the following result: every dense type of approx- 
imation appears on infinitely many discretizations of a generic homeomorphism. 
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Theorem 16. LetU he an open suhsef o/Homeo(X, A) andT he a type of approximation 
which is dense in U. Then for a generic homeomorphism f and for all Nq G there 
exists N > No such that fN ^ Tn- 

When we will try to obtain properties of discretizations of generic homeomorphisms, this 
theorem will be the second step of the proof, the first being to show that the approximation 
type we are interested in is dense. In section 6 we will try to establish as many corollaries of 
this theorem as possible, showing that the approach "prove that such approximation is dense 
then apply theorem 16" is a general method in the study of properties of discretizations^. 

To prove of theorem 16 we start from a dense type of approximations — in other words 
discrete applications — and we want to get properties of homeomorphisms. The tool that 
allows us to restore a homeomorphism from a finite map ajy * En E^-, i.e. do the 
opposite of a discretization, is the finite map extension proposition: 

Proposition 17 (Finite map extension). Let (xi, X2, . . . , x^) and • • • , ^fc) 'f^'^o 

k-tuples of pairwise distinct points in X \ dX (i.e. for all i ^ j, Xi ^ Xj and yi ^ yj). 
Then there exists f G Homeo(X, A) such that for all i, f{xi) = yi. Moreover f can he 
chosen such that if mdiXi d{xi^yi) < S then Id) < S. 

The idea of the proof is quite simple: it suffice to compose homeomorphisms whose 
support's size is smaller than e and which are central symmetries within this support. Then 
one constructs a sequence {zj)i<j<k such that zi = xi, yi = and d{zj^Zjj^i) < e/10. 
Composing k — 1 homeomorphisms as above, such that each one sends Zj on ^j+i, one 
obtains a conservative homeomorphism which sends xq on yo. Implementing these remarks 
is then essentially technical. A detailed proof can be found in section 2.2 of [Guil2]. 

With this proposition, we can build homeomorphisms from a dense type of approxima- 
tion: 

Lemma 18. Let U he an open suhset o/ IIomeo(X, A) and T = (T/v)ArGN he a type of 
approximation which is dense in U. Then for all Nq G the set of homeomorphisms f 
such that there exists N > Nq such that the discretization /at helongs to Tn 'is dense in U. 

Proof of lemma 18. Let f e Nq e N and e > 0. Since T is a type of approximation 
which is dense in there exists N > Nq and (Jat G Tn such that dN{f^cFN) < ^- Let 
Xi, . . . , Xq^ be the elements of En- Then for all ^, d{f{x£)^ cfn{xi)) < e. We modify aN into 
a bijection : En X whose discretization on En equals to aN' set cr^(xi) = cr(xi), 
cr^ is defined by induction by choosing <j^(x^) such that <j^(x£) is different from <j^(xi) 
for i < such that Pn{cf'n{^i)) = Pn{o'n{x£)) and that d{f{x£),a'^{x£)) < e. 

Since is a bijection proposition 17 can be applied to f{x£) and cr^(x£); this gives 
a measure-preserving homeomorphism (p such that (p{f{x£)) = cf'n{x() for all £ and such 
that (i((/:),Id) < e. Set = o /, we have d{f^f') < e and f'ixi) = cf'j^{xi) for all and 
therefore fN = o-n- D 

Proof of theorem 16. Let {xN,i)i<i<qN be the elements of En and consider the set (where 
Pat is the minimal distance between two distinct points oi En) 



This set is clearly a Gs set and its density follows directly from lemma 18. Moreover we 



^We will often take U = Homeo(X, A). 

^Although some properties can be proved in different ways: for instance corollary 37 can be proved 
in inserting horseshoes in the dynamics of a given homeomorphism, using the local modification theorem 
(theorem 30, for a presentation of the technique in another context see section 3.3 of [Guil2]); also a 
variation of corollary 25 can be shown in perturbing any given homeomorphism such that it has a periodic 
orbit whose distance to the grid E]\j (grids are assumed to refine) is smaller than the modulus of continuity 
of /, so that the actual periodic orbit and the discrete orbit fit together. 




ATqGN N>No 



can easily see that its elements satisfy the conclusions of the theorem. 



□ 
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5. Lax's theorem 

Now that we have shown theorem 16, we have to obtain dense types of approximation. 
Again, a theorem wih be systematicahy used: Lax's theorem, and more precisely its im- 
provement stated by S. Alpern [Alp76] which ahows us to approach any homeomorphism 
by a cychc permutation of the elements of a discretization grid. 

Theorem 19 (Lax, Alpern). Let f G Homeo(X, A) and e > 0. Then there exists Nq e N 
such that for all N > Nq, there exists a cyclic permutation cfn of En such that d{f^ (Jat) < e. 

As a compact metric set can be seen as a "finite set up to e", Lax's theorem allows us to 
see an homeomorphism as a "cyclic permutation up to e". In practice this theorem is used 
to "break" homeomorphisms into small pieces. Genericity of transitivity in Homeo(X, A) 
follows easily from Lax's theorem together with finite map extension proposition (again, see 
[APOO] or part 2.4 of [Guil2]). In our case, applying theorem 16, we deduce that infinitely 
many discretizations of a generic homeomorphism are cyclic permutations. The purpose of 
section 6 is to find variations of Lax's theorem (which are at the same time corollaries of 
it) concerning other properties of discretizations. 

We give briefiy the beautiful proof of Lax's theorem, which is essentially combinatorial 
and based on marriage lemma and on a lemma of approximation of permutations by cyclic 
permutations due to S. Alpern. Readers wishing to find proofs of these lemmas may consult 
section 2.1 of [Guil2]. This is here that we use the fact that discretizations are well 
distributed and well ordered. 

Lemma 20 (Marriage lemma). Let E and F he two finite sets and ^ a relation between 
elements of E and F. Suppose that the number of any subset E' d E is smaller than the 
number of elements in F that are associated with an element of E' , i.e.: 

yE' C E, Card(£;0 < Card {f e F \ 3e e E' : e ^ f} , 

then there exists an injection ^ : E ^ F such that for all e G E, e ^ ^{^)- 

Lemma 21 (Cyclic approximations in 6^, [Alp76]). Let q G N* and a e &q (&q is seen 
as the permutation group ofZ/qZ). Then there exists r e &q such that \r{k) — k\ < 2 for 
all k (where |.| is the distance in Z/qZ) and such that the permutation ra is cyclic. 

Proof of theorem 19. Let / G Homeo(X, A) and £ > 0. Consider A^o ^ N such that for 
all N > Nq, the diameter of the cubes of order N (given by the hypothesis "being well 
distributed") and their images by / is smaller than e. We define a relation ^ between cubes 
of order N — 1: C ^ C if and only if f{C) DC' ^ 9. Since / preserves A, the image of the 
union of i cubes intersects at least i cubes, so marriage lemma (lemma 20) applies: there 
exists an injection from the set of cubes of order N into itself (then a bijection) such 
that for all cube C, f{C) fl ^Ar(C) ^ 0. Let be the application that maps the center of 
any cube C to the center of the cube ^n{C), we obtain: 

dwifj ctn) < sup ( diam(C)) + sup ( diam(/(C))) < 2£. 

It remains to show that (Jat can be chosen as a cyclic permutation. We number the cubes 
such that two consecutive cubes are adjacent (this is possible because the grids are well 
ordered), then we use lemma 21 to obtain a cyclic permutation whose distance to ctat 
is smaller than 2 sup(diam((7)). Thus we have found Nq G 'N and for all N > Nq di cyclic 
permutation cr^ of En whose distance to / is smaller than Ae. □ 

6. Individual behaviour of discretizations 

Once we have set theorems 16 and 19, we can establish results concerning the behaviour 
of discretizations of a generic conservative homeomorphism. Here we study individual 
behaviour of discretizations, i.e. properties about only one order of discretization. As has 
already been said, applying theorem 16, it suffices to find dense types of approximation to 
obtain properties about discretizations. In practice, these dense types of approximations 
are obtained from variations of Lax's theorem (theorem 19). 
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Recall that the sequence (£^Ar)nGN of discretization grids is well distributed and well 
ordered (see definition 8), we denote by /n the discretization of an homeomorphism / and 
^(/at) the maximal invariant set of /at (i.e. the union of periodic orbits of /at). 

We will try to show that for simple dynamical properties (P) about finite maps and 
for a generic conservative homeomorphism /, infinitely many discretizations /at verify (P) 
as well as infinitely many discretizations verify its contrary. For instance, for a generic 
homeomorphism /, the maximal invariant set ^{/n) is sometimes as large as possible, 
i.e. ^{/n) = En (corollary 22), sometimes very small (corollary 25) and even better 
sometimes the number of elements of the image of En is small (corollary 37). In the same 
way stabilization time^^ is sometimes zero (corollary 22 for example), sometimes around 
Card(£;Ar) (corollary 27). Finally, concerning the dynamics of fN\Q(^fj^)^ sometimes it is a 
cyclic permutation (corollary 22) or a bicyclic permutation (remark 23, see also corollary 
35), sometimes it has many orbits (corollary 39). 



Firstly, we deduce directly from Lax's theorem that cyclic permutations of sets Ejy form 
a dense type of approximation in Homeo(X, A). Combining this with theorem 16, we obtain: 

Corollary 22 (Miernowski, [Mie05]). For a generic homeomorphism f G Homeo(X, A)^ 
for every Nq e'N, it exists N > Nq such that /at is a cyclic permutation}^. 

Remark 23. The same result is obtained for with bicyclic permutations, which are permu- 
tations having exactly two orbits whose lengths are relatively prime (see [Guil2]). 

The first variation of Lax's theorem concerns the approximation of applications which 
maximal invariant set is small. It is the contrary to what happens in corollary 22. 

Proposition 24 (First variation of Lax's theorem). Let f G Homeo(X, A). Then for all 
e^e' > 0, it exists A/'o G N such that for all N > Nq, it exists a map gn En En such 
that dN^f-, ctn) < £ cind 

Card(^](crAr)) _ Card(l](aAr)) 
Card(£^Ar) Qn 

and such that En is made of a unique (pre-periodic) orbit of gn- 

Proof of proposition 24- Let / G Homeo(X, A), e > and x a recurrent point of /. It 
exists r G N* such that d{x^f^{x)) < |; this inequality remains true for fine enough 
discretizations: it exists Xi G N such that if TV > Xi, then 

S ST 

d{x,XN)<-, d{f^{x)J^{xN))<- and — < . 
o o Qn 

Using the modulus of continuity of f^ and Lax's theorem (theorem 19), we obtain an 
integer Nq > Ni such that for all N > Nq, it exists a cyclic permutation gn of En such 
that dNif.cTN) < I and dN{f^,cr^) < §• Then 

d{xN,cr^N{^N)) < d{xN,x) ^d{xj^{x)) ^d{f^{x)J^{xN)) 
^d{f^{xN),cr^{xN)) 

£ 

We compose a at by the (non bijective) application mapping <j]y~^(xAr) on xn and being 
identity anywhere else (see figure 3), in other words we consider the application 

^/ ^\xn iix = cr]^^{xN) 

^ \crN{x) otherwise. 

The map has a unique injective orbit whose associated periodic orbit ^^(cr^) has length r 
(it is (xat, crAr(xAr), . . . , cf^^^^xn)))- Since d{f, cr'j^) < e, the map cr^ verifies the conclusion 
of the proposition. □ 
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Figure 3. Modification of a cyclic permutation 



A direct application of theorem 16 gives: 
Corollary 25. For a generic homeomorphism f G Homeo(X, A)^ 

Card(l^(/^)) 
Card(E^) 

Specifically for all £ > and all Nq G exists N > Nq such that ^cardiEN)^ ^ ^ 

5iic/i ^/la^ £^Ar is made of a unique (pre-periodic) orbit of f^- 

The same kind of idea leads to the following proposition: 

Proposition 26. Let f G Homeo(X, A) having at least one periodic point of period p. 
Then for all e > 0, it exists Nq e N such that for all N > Nq, it exists an application 
gn ' Ejsf Ejsf with Card(l](crAr)) = p and dN{f^crN) < s, such that En is made of a 
unique (pre-periodic) orbit of cfn- 

Proof of proposition 26. Simply replace the recurrent point by a periodic point of period p 
in the proof of proposition 24. □ 

We will show that owning a periodic point is a generic property in Homeo(X, A) (see also 
[DFOO] or part 3.2 of [Guil2]), a proper application of Baire's theorem leads to a refinement 
of corollary 25: 

Corollary 27. Let e > 0. A generic homeomorphism f G Homeo(X, A) has a periodic 
point whose orbit is e-dense. Moreover, f^ has a unique periodic orbit for infinitely many 
integers N , whose period does not depend on N e 'N and equals the smallest period of 
periodic points of f which are e-dense; moreover En is covered by a single (pre-periodic) 
orbit of fN- 

To prove corollary 27 we have to introduce the concept of persistent point. 

Definition 28. Let / G E[omeo(X). A periodic point x of f with period p is said persistent 
if for all neighbourhood U of x, there exists a neighbourhood V of / in Homeo(X) such 
that every / G V has a periodic point x e U with period p. 



^^l.e. the smallest integer k such that /^(En) = Q(fN). 

^^In fact, T. Miernowski proves "permutation" but his arguments, combined with the version of the 
Lax's theorem we gave, show "cyclic permutation". 
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Example 29. The endomorphism h = Diag(Ai, . . . , A^) of R^, with fl = ^ K 7^ ^ for 
ah is measure-preserving and has a persistent fixed point at the origin (see e.g. [KH95], 
p. 319). Let 5 be a reflection of R^, the application h o f is also measure-preserving and 
has a persistent fixed point at the origin. 

Finally, we use the theorem of local modification of conservative homeomorphisms, which 
allows us to replace locally a homeomorphism by another. Although it is geometrically "ob- 
vious" and it has an elementary proof in dimension two, it is easily deduced from (difficult) 
annulus theorem. For more details one may refer to [DFOO] or part 3.1 of [Guil2]. 

Theorem 30 (Local modification). Let cfi, (72, ti and T2 he four bicollared embeddings^^ 
o/S^-^ in such that <Ji is in the bounded connected component^^ of cf2 and ri is in 
the bounded connected component of T2- Let Ai be the bounded connected component of 
R"' \ (Ji and Bi the bounded connected component of R^ \ ri ; H be the connected component 
o/R"' \ (cTi LJcr2) having aiUa2 as boundary and A the connected component o/R"' \ (ri Ur2) 
having ri U r2 as boundary; A2 be the unbounded connected component of R^ \ <J2 and B2 
the unbounded connected component ofW^ \ T2 (see figure 4)- 

Suppose that Leb(74i) = Leh{Bi) and Leb(S) = Leb(A). Let fi : Ai ^ Bi be two 
measure-preserving homeomorphisms such that either each one preserves the orientation, 
or each one reverses it. Then it exists a measure-preserving homeomorphism f : R"' R^ 
whose restriction to A\ equals fi and restriction to A2 equals f2- 

Since this theorem is local, it can be applied to an open space O instead of R^, or even 
better, together with Oxtoby-Ulam theorem (theorem 12), to any domain of chart of a 
manifold X instead of R"^ and measure A instead of Lebesgue measure 

Proof of corollary 27. Let e > and Up the set of homeomorphisms that have at least one 
periodic persistent point with period p whose orbit is £-dense (strictly), but no periodic 
point with period less than p whose orbit is e-dense (not strictly). The sets Up are open 
and pairwise disjoints. Moreover their union is dense. Indeed, for / G Homeo(X, A), if / 
has at least one periodic point whose orbit is e-dense, then we choose a periodic point x 
whose orbit is e-dense and whose period p is minimal. Then we perturb / such that this 
point is persistent, avoiding the creation of periodic points of smaller periods: we choose a 
neighbourhood D of x such that the sets f{D), • • • , f^~^{D) are pairwise disjoints and 
we replace locally / by the map ho f~^P~^^ (where h is one of the two maps of the example 
29, depending of whether /~(^~^) preserves orientation or not) in the neighbourhood of 
/~^(x), using the theorem of local modification (theorem 30), such that / do not change 
outside the union of the sets f^{D). On the contrary, if / does not have any periodic point, 
we pick a recurrent point and perturb / such that this point is periodic, and then we apply 
the previous case. 

"'^^An embedding i of a manifold M in is said bicollared if it exists an embedding j : M x [— 1, 1] — ^ 
such that Jmx{o} — 

""^^By Jordan-Brouwer theorem the complement of a set which is homeomorphic to S^"-*^ has exactly 
two connected components: one bounded and one unbounded. 
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Let p G N*. Proposition 26 indicates that the set of maps whose longest periodic orbit 
has length p is a dense type of approximation in Up. So, by theorem 16, it exists open sets 
Op^N such that Hatgn ^p,n is a dense Gs of Up, made of homeomorphisms of Up whose an 
infinity of discretizations have only one injective orbit whose associated cycle has length p 
and is e-dense. 

Set On = UpGN* ^p,N- Since the sets Up are pairwise disjoints. 

The right side of the equation forms a dense subset of Homeo(X, A), so the left side is a 
dense Gs of Homeo(X, A). It is made of homeomorphisms / such that for infinitely many 
A/" G N, the grid En is covered by a single forward orbit of /at whose associated periodic 
orbit is e-dense and whose length is the shortest period of periodic points of / whose orbit 
are e-dense. □ 

Remark 31. The shortest period of periodic points of generic homeomorphisms has not 
any global upper bound in Homeo(X, A): for example, for all p G N, there is a open set 
of homeomorphisms of the torus without periodic point of period less than p (e.g. the 
neighbourhood of a nontrivial rotation) and this property remains true for discretizations. 

We now want to get a discrete analogue to topological weak mixing (which is generic in 
Homeo(X, A), see e.g. [APOO] or part 2.4 of [Guil2]). 

Definition 32. A homeomorphism / is said topologically weakly mixing if for all nonempty 
open sets {Ui)i<M and {U[)i<Mi it exists m G N such that f^{Ui) fl U!- is nonempty for all 
i < M. 

The proof of the genericity of topological weak mixing starts by an approximation of 
every conservative homeomorphism by another having £-dense periodic orbits whose lengths 
are relatively prime. The end of the proof lies primarily in the use of Baire's theorem 
and Bezout's identity. In the discrete case, the notion of weak mixing is replaced by the 
following: 

Definition 33. Let £ > 0. A finite map a at is said e -topologically weakly mixing if for all 
M G N and all balls {Bi)i^M and {B'-)i<^M with diameter £, it exists m G N such that for 
all i 

a^{B,nEN)n{B'^nEN)^9. 

The first step of the proof is replaced by the following variation of Lax's theorem: 

Proposition 34 (Second variation of Lax's theorem). Let f G Homeo(X, A) be a homeo- 
morphism whose all iterates are topologically transitive. Then for all e > and all M G N*^ 
there exists TVq G N such that for all N > Nq, it exists gn ' En En which has 
M s-denses periodic orbits whose lengths are pairwise relatively prime, and such that 
dNif.cTN) < e. 

Proof of proposition 34- We prove the proposition in the case where M = 2, the other 
cases being easily obtained by an induction. Let s > and / be an homeomorphism whose 
all iterates are topologically transitive. Then it exists G X and p G N* such that 
{xo, . . . ^ f^~^ (xq)} is e-dense and d{xQ, f^{xo)) < e/2. Since transitive points of form 
a dense Gs subset of X, while the orbit of xq form a F^r set with empty interior, the set 
of points whose orbit under is dense and disjoint from that of xq is dense. So we can 
pick such a transitive point y^. Set yi = f{yo)' Then there exists a multiple qi of p such 
that the orbit {yi^ - • - ^ f^^~^{yi)} is e-dense and d{yi, f^^{yi)) < e/2. Again, by density, 
we can chose a transitive point y2 whose orbit is disjoint from that of xq and ^i, with 
d{yi,y2) < s/2 and d{yo,yi) — d{yo,y2) > e/4. Then it exists a multiple q2 of p such that 
<^(^2, /^^(^2)) < s/2. And so on, we construct a sequence {ym)i<m<£ such that (see figure 
5): 

(i) for all m, it exists Qm > such that p\qm and d{ym^ f^'^iym)) < ^/2, 
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Figure 5. Construction of the sequence {ym)i<m<£ 

(ii) the orbits {xq, . . . , /^~H^o)} and {ym, • • • , /^"""H^m)} {m going from to ^ - 1) are 
pairwise disjoints, 

(iii) for ah m, d{ym,ym+i) < e/2 and d{yo,ym) - d{yQ,ym+i) > ej^, 

(iv) yi = yQ. 

Let (Jat be a finite map given by Lax's theorem. For ah N large enough, satis- 
fies the same properties (i) to (iii) than /. Changing gn at the points (J^~^{{ym)N) 
and (j^^ {{xo)n)^ we obtain a finite map such that cr'js^^^ {{ym)N) = {ym-\-i)N and 
cr%{{xo)N) = {xo)n. Thus the orbit of {xo)n under is 2£:-dense and has period p 
and the orbit of {yo)N under cr^ is 2£-dense, disjoint from which of {xo)n and has period 
1 + ^1 + • • • + Q£-i relatively prime to p. □ 

Corollary 35. For a generic homeomorphism f G Homeo(X, A)^ for all s > and all 

Nq e it exists N > Nq such that /n is e-topologically weakly mixing. 

Proof of corollary 35. Again, we prove the corollary in the case where M = 2, other cases 
being easily obtained by induction. Let £ > and Nq G N. All iterates of a generic 
homeomorphism / are topologically transitive: it is an easy consequence of the genericity 
of transitivity (see e.g. corollary 22 or theorem 2.11 of [Guil2]); we pick such a homeo- 
morphism. Combining theorem 16 and proposition 34, we obtain N > Nq such that /at 
has two e/3-dense periodic orbits whose lengths p et q are coprime. We now have to prove 
that /at is e-topologically weakly mixing. Let 5i, B[ and B2 be balls with diameter 
£. Since each one of these orbits is e/S-dense, it exists xat G X which is in the intersection 
of the orbit whose length is p and 5i, and yN ^ X which is in the intersection of the orbit 
whose length is q and B2. Similarly, it exists two integers a and h such that f^ixjy) G B[ 
and /^(^at) e B'^. 

Recall that we want to find a power of /at which sends both x at in 5^ and y^ in B2. It 
suffices to pick m G N such that m = a-\-ap = b-\- j3q. Bezout's identity states that there 
exists two integers a and {3 such ap — /3q = b — a. Set m = a-\- ap, adding a multiple of pq 
if necessary, we can suppose that m is positive. Thus f^{xN) G B[ and f^{yN) ^ ^2- ^ 

We now establish a new variation of Lax's theorem. It asserts that every homeomorphism 
can be approximated by a finite map whose image has a small cardinal, unlike what happens 
in Lax's theorem. 

Proposition 36 (Third variation of Lax's theorem). Let f G Homeo(X, A) and 1^ : N ^ 
a map which tends to +00 in +00. Then for all e > 0, it exists A/'o G N such that 
for all N > Nq, it exists a map - such that Card((TAr(£^Ar)) < ^{N) and 

dNif.cTN) < e. 

Proof of proposition 36. Let / G IIomeo(X, A), 'i^ : N ^ a map which tends to +00 
in +00 and £ > 0. By Lax's theorem (theorem 19) it exists A"! G N such that for all 
A^ ^ A'l, it exists a cyclic permutation cfn : E^ En whose distance to / is smaller than 
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For N > A^i, set = P/v^ o a^- Increasing A^i if necessary we have d{f^a^j^) < 
regardless of N. Moreover Card(cr^(£^Ar)) < Qni, if we choose A^o large enough such that 
for ah TV > TVo we have qn^ < ^(7V), then Card(cr^(£;Ar)) < i^{N). We have shown that 
the map satisfies the conclusions of proposition for all N > Nq. □ 

Corollary 37. Let : N ^ a map which tends to +00 in +00. Then for a generic 
homeomorphism f G Homeo(X, A)^ 

hm Card(/.(^^))^o. 
In particular, generically, lim ^^^^^^.'^^^^y^^ = 0. 

Proof of corollary 37. Remark that if we replace d^N) by ^/l9{N)^ it suffices to prove that 
for a generic homeomorphism, lim ^^^^^J^^^^^^ < 1. It is easily obtained in combining 
theorem 16 and proposition 36. □ 

So far all variations of Lax's theorem built finite maps with a small number of orbits. 
With the additional assumption that the sequence of grids is self-similar, we show a last 
variation of Lax's theorem approaching every homeomorphism with a finite map with a 
large number of orbits. 

Proposition 38 (Fourth variation of Lax's theorem). Assume that the sequence of grids 
{EN)Ne'N 'is self-similar. Let f G Homeo(X, A) and 'i^ : N ^ R such that 'd{N) = o{qN)- 
Then for all e > it exists A^i G N such that for all N > Ni, it exists a permutation 
of En such that d^if^crN) < s and that the number of cycles of is greater than 
i9{N). Moreover these i9{N) cycles of aN are conjugated to a cyclic permutation of Enq by 
bisections whose distance to identity is smaller than e. 

Proof of proposition 38. Let e > 0, for all A/'o G N large enough, Lax's theorem gives us a 
cyclic permutation of E^q whose distance to / is smaller than e. Since the grids are 
self-similar, it exists A^i G N such that for all N > Ni, the set En contains Qn/qnq ^ i^{N) 
disjoint subsets each one conjugated to a grid Enq by a bijection hj whose distance 
to identity is smaller than e. On each Ej^, we define a at as the conjugation of cr^^ by hj; 
outside these sets we just pick aN such that dN^f^cFN) < ^- Since the distance between hj 
and identity is smaller than we have dN{f^cFN) < 26. Moreover, aN has at least i9{N) 
cycles; this completes the proof. □ 

The application of theorem 16 gives us: 

Corollary 39. We still assume that the sequence of grids {EN)Ne'N is self-similar. Let 
i9 '.N such that i9{N) = o{qN)- Then for a generic homeomorphism f G Homeo(X, A) 
and for infinitely many integers N, the discretization fN of f have at least i^{N) cycles 
which are pairwise conjugated. 



7. Average behaviour of discretizations 

We now want to study the average behaviour of discretizations of a generic homeomor- 
phism. For example one could imagine that even if for a generic homeomorphism /, the 
event "/at is a cyclic permutation" appears for infinitely many orders, it is still quite rare. 
More precisely, we study the frequency of occurrence of properties related to the discretiza- 
tions of generic homeomorphisms in the Cesaro sense: given a property (P) concerning 
discretizations, what is the behaviour of the proportion between 1 and M of discretizations 
satisfying the property (P), when M goes to infinity? For this study, we assume that the 
sequence of discretization grids refines (which is true for example for discretizations upon 
uniform grids of orders powers of an integer, see section 3). This prevents us from tricky 
arithmetic problems about overlay of grids. 
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Definition 40. Let / G Homeo(X, A). We say that a property (P) about discretizations 
is satisfied in average if for all Nq e N and all £ > 0, it exists N > Nq such that the 
proportion of integers M G {0, . . . , N} such that /m satisfies (P) is greater than 1 — i.e. 

lim -J— Card |M G {0, . . . , A/"} I /m satisfies (P)| = 1. 

We will show that most of the dynamical properties studied in the previous section 
are actually satisfied on average for generic homeomorphisms. To start with we set out a 
technical lemma: 

Lemma 41. Let T be a dense type of approximation in Homeo(X, A). Then for a generic 
homeomorphism f G Homeo(X, A)^ for all e > and all a > 0^ the property (P) .* "En 
contains at least a disjoints subsets which fulfills a proportion greater than 1 — e of En, 
each one stabilized by /at and such that the restriction of f^ to each one is conjugated to a 
map of T by a bijection whose distance to identity is smaller than e^^ is satisfied in average. 

In practice, this lemma provides many properties satisfied on average, for instance: 

• quantitative properties on discretizations, such as owning at least M periodic orbits, 

• properties of existence of sub-dynamics on discretizations, such as owning at least 
one dense periodic orbit. 

Proof of lemma J^l. Let us consider the set 

C _ r\ \ \ \ f ^ Homeo(X, A) | 1 

Li N>N \N^iC^'^{Me{0,...,N}\fMS^tis^es{P)}>l-ej' 

We want to show that C contains a dense Gs of Homeo(X, A). The set C is a of the 
generic set Hatgn ^at, it suffices to prove that it is dense in Vn- Let / G Homeo(X, A), 
Xo G N, (5 > and e > 0. To prove it we want to find a homeomorphism g whose distance 
to / is smaller than S and an integer N > Nq such that 

Card {M G {0, . . . , X} I satisfies (P)} > 1 - £. 

It is simply obtained in combining the density of the type of approximation T and the fact 
that the grids refines. □ 

This corollary allows us to obtain some properties about the average behaviour of dis- 
cretizations. For instance there is an improvement of corollary 27: 

Corollary 42. For a generic homeomorphism f G Homeo(X, A)^ the property "f^ has a 
e-dense periodic orbit and the cardinal ofQ{fN) is smaller than 'd{N) = o^qn)^^ is satisfied 
in average. 

Or an improvement of corollary 35: 

Corollary 43. For a generic homeomorphism f G IIomeo(X, A) and for all e > {), the 

property "f^f is e-topologically weakly mixing (see definition 33)^^ is satisfied in average. 

Or even an improvement of corollary 37: 
Corollary 44. For a generic homeomorphism f G Homeo(X, A) and for all e > 0^ the 
property ^^ ^^cardCg^T^^ < s^^ is satisfied in average. 

And an improvement of corollary 39: 

Corollary 45. For a generic homeomorphism f G Homeo(X, A) and for all M G 

property '^f^ has at least M periodic orbits^' is satisfied in average. 

However, note that the most simple property about discretizations, i.e. being a cyclic 
permutation, can not be proved by using lemma 41. To do this, we need a slightly more 
precise result, that we will not prove here. 

Proposition 46. For a generic homeomorphism f G Homeo(X, A); the property '^f^ is a 
cyclic permutation^^ is satisfied in average. 



20 



PIERRE-ANTOINE GUIHENEUF 



Remark 47. However, the property of approximation by bicyclic permutations in average 
can not be proven with this technique. 

Remark 48. A simple calculation shows that everything that has been done in this section 
also applies to the behaviour of discretizations in average of Cesaro average, in average of 
average of Cesaro average etc., i.e. when studying quantities 

JTTl i JTTl ^^'^ ^ {0- • - ^1} I /m satisfies (P)}, 

]^ E ]^ E ^Card {M e{0,...,iV,} I /m satisfies (P)}... 

N2=0 ^ Ari=o ^ 

8. Behaviour of all the discretizations 

In the previous sections we showed that the dynamical behaviour of discretizations 
seems very chaotic depending of the order of discretization. In contrast, the dynamics 
of a generic homeomorphism is well known (see e.g. [Guil2]) and independent from the 
homeomorphism. One even has a 0-1 law on Homeo(X, A) (see [GK98] or the last chapter of 
[Guil2]) which states that either a given ergodic property on conservative homeomorphisms 
is generic, or its contrary is generic. Thus, the dynamics of a generic homeomorphism and 
that of its discretizations seem quite of independent. In fact one can deduce some dynamical 
features of a generic homeomorphism from the corresponding dynamical features of all its 
discretizations. For instance, the following property can be easily deduced from corollary 
27. 

Proposition 49. Let f be a generic homeomorphism o/Homeo(X, A) and p be an integer. 
Then f has a periodic orbit with period p if and only if there exists infinitely many integers 
N such that f^ has a periodic orbit with period p. 

Proof of proposition 49. An easy variation of corollary 27 shows that if p is a period of a 
periodic orbit of a generic homeomorphism /, then there exists infinitely many discretiza- 
tions /at such that /at has a periodic orbit with period p. The other implication of the 
proposition follows easily from a compactness argument. □ 

We now try to obtain information about invariant measures of a generic homeomorphism 
from invariant measures of its discretizations. More precisely, given all the invariant mea- 
sures of discretizations of a generic homeomorphism, what can be deduced about invariant 
measures of the initial homeomorphism? A first step in this study was performed by T. 
Miernowski in 2006 in part 8 of his article [Mie06] : 

Proposition 50 (Miernowski). Let f : X ^ X be a uniquely ergodic homeomorphism and 
ji^ his unique invariant probability measure. For all N e N let C En be a periodic 
cycle of fN cind the uniform probability measure on jn- Then the weak convergence 
1^^ occurs independently of the choice of the cycles • 

The proof of this proposition essentially consists in an appropriate application of Prokho- 
rov's theorem which express the compactness of the set of probability measures on X. We 
now set a theorem of this kind for generic homeomorphisms. Recall that /i^^ is the limit in 
the sense of Cesaro of the pushforwards by iterates of fjy of uniform measure on En (see 
definition 6): 

^ m—l 

I^N= 1™ — y](/Ar)lAAr. 
m^oo rn ^ — ^ 

The measure is supported by the maximal invariant set of /at; it is uniform on every 
periodic orbit and the total weight of a periodic orbit is proportional to the size of its 
basin of attraction. The following theorem expresses that we can obtain all the invariant 
measures of a generic homeomorphism from all the invariant measures of its discretizations: 
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Theorem 51. Let f G Homeo(X, A) be a generic homeomorphism and suppose that the 
sequence of grids {EN)Ne'N is self-similar. Let Mn the set of probability measures on 
En that are invariant under fN- Then the upper limit over N (for Hausdorff metric) of 
the sets M.n is exactly the set M. of probability measures that are invariant under f . 

Before giving a detailed proof of theorem 51 let us give its main arguments. To show that 
the upper limit of sets A^at is a subset oi by a compactness argument, we only have 
to prove that limit points of a sequence (z/^)ArGN, such that every measure is invariant 
under /at, are invariant under /. This can be easily verified in using uniform convergence of 
/at to / and equicontinuity of the measures v^^. Then (again by a compactness argument) 
we have to show that generically every invariant measure is a limit point of a sequence of 
discretizations. An ad hoc application of Baire's theorem reduces the proof to which of the 
following variation of Lax's theorem^^: 

Lemma 52 (Ergodic variation of Lax's theorem). Suppose that the sequence of grids 
{EN)Ne'N 'is self-similar. For all f G Homeo(X, A)^ for all collection of f -invariant mea- 
sures . . . for all € > and /cq, Nq G N.* 

(1) It exists g G Homeo(X, A) and N > Nq such that d{f^g) < e, (i(i^i,/i^) < ^ and 
fi^j^ is g -invariant. 

(2) It exists g G Homeo(X, A) and N > Nq such that d{f^g) < e, and for all i < £ it 
exists a measure ^ such that d{vi^vf j^) < ^ and that is invariant under g 
and gN- 

Suppose first that i = 1 and that ui is ergodic. For this purpose we apply Birkhoff 's 
theorem to f^ui^cp and a recurrent point x: for all M large enough, 

1 ^"^ r 

k=0 

Since x is recurrent we can choose an integer M large enough such that x f^ (x). First 
we approximate / by a cyclic permutation ctat given by Lax's theorem, then we slightly 
modify into a map by choosing cr^(crjjf (xat)) = xn, as in proposition 24. The 

measure is the uniform measure on the orbit xat, . . . , ajjf "^(xat), so it is close to ui 
with respect to the test function (p. This property can be obtained for all test functions by 
a diagonal-like method. The proof of the lemma when ui is not ergodic but only invariant is 
obtained by approximating the invariant measure by a finite convex combination of ergodic 
measures; this proves the first point of the lemma. To obtain the latter, it suffices to use the 
self-similarity hypothesis: on at least i sub-grids we apply the first point to each measure 

Proof of theorem 51. For all 5 > ^, for all N large enough, Mn ^ B{M^S). For a set A 
and J > 0, we note B{A^S) the set of elements whose distance to A is smaller than S. This 
inclusion follows easily from the upper semi-continuity of the application g ^ M.^ (where 
M.^ denotes the set of Borel probability measures which are invariant under g) and the 
compactness of V. 

Generically, for all S > and all Nq G 'N, it exists N > Nq such that M C B{Mn^ 
Recall that we denote by V the set of all Borel probability measures on X. Set d a distance 
on V defining the weak-* topology. Thereafter homeomorphisms will be taken in the generic 
set 

n 

ATGN 

made of homeomorphisms whose Nih discretization is uniquely defined for all N . Consider 

(Aro,feo)GN2 



The first point of this lemma is somehow a weak version of the latter. It will be useful in the next 
section. 
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where 



I iVeN 



3N > No s.t. M C B{Mn, |-) 



Trivially, if / G then the upper limit of the sets Mn contains M. 

To show that the sets ONo,ko open, it suffices to remark that if M C B(^Mn^ ^) for 
a given then it is also true on a neighbourhood of / since g M% is constant on a 
neighbourhood of / and g ^ M.^ is upper semi-continuous. 



, TO /-» V"! C 1 -1- r\ 



It remains to show that the sets OnqM are dense; it follows from the second point of 
the variation of Lax's theorem (lemma 52): by upper semi-continuity of ^ M.^ ^ one 
has C B(^M^ ^ ^) for all g close enough to /. Moreover, by compactness, it exists 
/-invariant measures vi^. . . such that (Z{j^B{yi^^Y So it suffices to find g close 
to / such that for a > A^o, it exists ^7v-invariant measures . . . such that 

d{vi^vf < ^ for all i, but this is exactly the conclusion of the lemma. □ 

Proof of lemma 52. To begin with we prove the first point of the lemma. Suppose first 
that the measure v = vi \s ergodic. We want to show that there exists a homeomorphism 
g whose distance to / is smaller than e, and an integer N > Nq such that d{iy^ jj^^jy) < 
Let {(pj)j^-[si be a sequence which is dense in the set of continuous maps from X into R. 
By Prokhorov's theorem, there exists i G N such that if we have 

1 

< - 
i 



(fj du — (fj d/i 



for all j < i, then d{h'^ /i) < ^. 

Since u is ergodic, for all continuous map (f, by Birkhoff 's theorem. 



M-l 
m=0 



for z/-a.e. x. Let x G X be a recurrent point for / satisfying equation (1) for all (fj (such 
points form z/-full measure set). For all j there exists Mj G N such that for all M > Mj 
one has 

M-l 

- ^^,o/-(x)- /^,d 



Let M- = maxj<i Mj. Since x is recurrent, there exists r > M- such that d{x, /^(^)) ^ 
e/4. Let (Jat be a map from into itself given by Lax's theorem: it is a cyclic permutation 
and its distance to / is smaller than e/2. For all TV large enough the orbit (cr]^(xAr))o<m<r 
shadows the orbit {f^{x))o<rn<T^ thus d{xN^o-]^{xN)) < e/2. Then we "close" the orbit of 
xn between the points xn and o-J^{xn), i.e. we set (as in proposition 24, see also figure 3) 



xn if Vn cr]v ^{^n) 

o'Nil/N) otherwise. 



Then dN{a'j^^f) < e and (xat, cr^~^(^Ar)) is a periodic orbit for whose basin of 
attraction is the whole set E^. 

Since the periodic orbit (cr^(x))o<m<r attracts En^ for all M' large enough and yN G 
En we have 

— ^ ^, o a'^iy^) ^-^^ - E ° "nM- (3) 

m=0 m=0 



With proposition 17, the same way as in lemma 18, we construct a homeomorphism g 
from the map such that the discretization {g^)N and fit together, and such that 
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and the restriction of {g )n to the orbit of xn fit together. Since (iAr(/, cr^) < we can 
furthermore assume that d{f^g^) < e. Thus, equation (3) imphes that 

r-1 



1 



(4) 



m=0 



Since ipj is continuous, increasing N (if necessary) such that the orbit (crX^(xAr))o<m<r 
shadows (/^(a:))o<m<r better, we have for ah j < i: 



m=0 



m=0 



< 



2i 



(5) 



We now have ah the necessary estimations to compute the distance between /i^ and u: 
by equation (4) and triangle inequahty. 



r-1 



/ (fj d/j.% - f (fjdu = - ^ (pj o a^{xN) - I (Pj du 

1 ^"^ 1 ^"^ 



=0 
r-1 



Hence, with equation (2) apphed to M = r and equation (5), 

1 



d/i 



N 



Lpj dv 



< 



We obtain , i^) < Since /i^ is invariant under (g^ is a cychc permutation of 

the orbit of x at), is the desired homeomorphism. 

In the general case the measure v is only invariant (and not ergodic). It reduces to 
the ergodic case by the fact that the set of invariant measures is a compact convex set 
whose extremal points are exactly the ergodic measures: by Krein-Milman theorem, for 
all M > 1 there exists an /-invariant measure which is a finite convex combination of 
ergodic measures: 

and whose distance to v is smaller than 7^. 

Then we use the hypothesis of self-similarity of the grids: for all e > and all P G N 
there exists a map ctn whose distance to / is smaller than e and whose at least p > P 
cycles are conjugated to a permutation of a coarser subdivision, via the composition of a 
fixed homothety by a translation. Moreover the union on these p cycles fulfills the set En 



with a proportion of at least 1 
and r we set 



and 



2fe, 



^. Take P G N such that ^ < 



2ko ' 



For all j between 1 



a: 



A[) = l 



1 we can 



Then for all j > we have \Xj — | < and |Ao — AqI < Since XI A^- 
associate bijectively a set Aj made of pXj cycles to each integer j; the union of the sets 
Aj fills the set En with a proportion of at least 1 — 2^^- We do the previous construction 
concerning ergodic measures for each j: increasing N if necessary, the cycles of each Aj 
can be changed such that each of them carries a (unique) invariant probability measure 
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satisfying d{iyj,h'j') < These modifications transform the map into a map 
satisfying d{(7'^^(7N) < Then the measure ji^^ associated with the map cr'^ satisfies: 



and we have: 



2fco 

3=0 



j=0 j=0 ° 



j=0 j=0 

2r 1 1 3 

< \ \ < . 

- P 2ko 2ko ~ 2ko 

Thus d{v^ii^^) < As above the desired homeomorphism is obtained by extending 
to a homeomorphism of X. We have proven the first point of the lemma. 

Finahy we handle the second point of the lemma. To do that, it suffices to apply the 
first point on I different sub-grids of a given grid Enq- On the i-th sub-grid, define the 
discretization such that it has a single invariant measure which is close to vi. This 
proves the lemma. □ 

There is a kind of duality between invariant measures and invariant compact sets, thus 
the previous theorem is also true for invariant compact sets: 

Proposition 53. Let f G Homeo(X, A) be a generic homeomorphism and suppose that 
the sequence of grids {EN)Ne'N 'is self-similar. Set Xn the set of subsets of that are 
invariant under /at. Then the upper limit over N (for Hausdorff metric) of the sets Xn is 
exactly the set of compacts that are invariant under f . 

Sketch of proof of proposition 53. To show that all limit point iiT of a sequence {K^) of 
compact invariant sets is a compact invariant set, it suffices to calculate dniK^ fi^))- For 
the other inclusion, as in the case of invariant measures, it follows from an approximation 
lemma whose proof is the same as for invariant measures: □ 

Lemma 54 (Compact variation of Lax's theorem). Suppose that the sequence of grids 
{En)n^-^ is self-similar. For all f G Homeo(X, A)^ for all collection of f -invariant com- 
pacts i^i, . . . , K£, for all e > and ko^ Nq G N.* 

(1) It exists g G Homeo(X, A) and N > Nq such that d{f^g) < e, <i//(i^i, 1^(/Ar)) < ^ 
and ^{fN) 'is g -invariant. 

(2) It exists g G Homeo(X, A) and N > Nq such that d{f^g) < e, and for all i < i, 
it exists a compact Kf such that dniKi^Kf j^) < CLTid that Kf ^ is invariant 
under g and gN- 

9. Physical measures 

In the previous part we proved that generically, the upper limit of the sets of invariant 
measures of discretizations is the set of invariant measures of the initial homeomorphism. 
It expresses that the sets of invariant measures of discretizations behave "as chaotically as 
possible". However, one might expect that physical measures (Borel measures fi such that 
/i = /i^ is verified on a set of points x with A positive measure, see e.g. [You02]) play a 
specific part: their definition expresses that they are the measures that can be observed in 
practice, because they governs the ergodic behaviour of A a.e. point. So one can hope that 
the natural invariant measures /ij^ of /at, which can be seen as the physical measures of 
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/at, converge to the physical measures of /. This expectation is supported by the fohowing 
variation of proposition 50, obtained by replacing by if f is uniquely ergodic, then 
the measures ji^^^ converge weakly to the only measure ji^ that is invariant under f. 

In this section we show that this is not at all the case: the sequence of measures (/i^^) 
accumulates of the whole set of /-invariant measures^^. More precisely one has the following 
theorem: 

Theorem 55. If the sequence of grids {EN)Ne'N is self-similar, for a generic homeomor- 
phism f G Homeo(X, A)^ the set of limit points of the sequence (/i-j^)]^^-^ is exactly the set 
of f -invariant measures. 

This theorem can be seen as a discretized version of the following conjecture: 

Conjecture 56 (F. Abdenur, M. Andersson, [AA12]). A homeomorphism f which is 
generic in the set of homeomorphisms of X (without measure preserving hypothesis) that do 
not have any open trapping set is wicked^ i.e. it is not uniquely ergodic and the measures 



accumulate on the whole set of invariant measures under f. 

The behaviour described in this conjecture is the opposite of that consisting of possessing 
a physical measure. 

Sketch of proof of theorem 55. The proof is similar to which of theorem 51: the set A is 
replaced by 



Remark 57. Taking over the proof of the theorem, if we set x ^ then generically 

the measures ^ accumulate on the whole set of /-invariant measures. This seems to con- 
tradict the empirical observations made by A. Boyarsky in 1986 (see [Boy86] or [GB88]): 
when a homeomorphism / has a unique ergodic measure /i-^ which is absolutely continu- 
ous with respect to Lebesgue measure "most of" the measures tend to measure ji^ . 
However, the author does not specify in what sense he means "most of the points", or if his 
remark is based on a tacit assumption of regularity for /. 

Note that as in the previous section, we have a compact counterpart of theorem 55: 

Proposition 58. // the sequence of grids {EN)Ne'N is self-similar, for a generic homeo- 
morphism f e Homeo(X, A)^ the set of limit points of the sequence ((7(/iv))ArGN is exactly 
the set of f -invariant compact sets. 

Part 2. Discretizations of a generic dissipative homeomorphism 

Throughout this section, we fix a manifold X and a good measure A as defined in section 
2. Now homeomorphisms are no longer supposed conservative. More formally, we recall that 
we note Homeo(X) the set of all dissipative homeomorphisms of X (without assumption 
of conservation of a given measure), we are interested in properties of discretizations of 
generic elements of Homeo(X). The manifold X is equipped with a sequence {En)n^-^ of 
discretization grids. 




m— 1 




and the second point of lemma 52 is replaced by the first. 



□ 



It is easy to see that if {i'N)Ne'N ^ sequence of measures, with ujy invariant under fjsj for all N, then 
its limit set is included in the set of /-invariant measures. Again, this reflects the fact that the behaviour 
of discretizations of a generic homeomorphism is "as chaotic as possible". 
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Generic topological properties of dissipative homeomorphisms have been discussed in 
a survey of E. Akin, M. Hurley and J. Kennedy [AHK03]. However, ergodic properties 
of generic homeomorphisms had not been studied until the work of F. Abdenur and M. 
Andersson in [AA12], in which the authors are interested in the typical behaviour of Birkhoff 
averages. For this purpose they establish a technical lemma they call shredding lemma 
which allows them to show that a generic homeomorphism is weird. 

Definition 59. A homeomorphism / is said weird if almost every point x G X (for A) has 
a Birkhoff 's limit /i^, and if / is totally singular {i.e. there exists a Borel set with total 
measure whose image by / is null measure) and does not admit any physical measure. 

In opposition to the conservative case, the lack of conservation of a given measure allows 
us to perturb any diffeomorphism to create attractors whose basins of attraction cover 
almost the entire space X. This is an open property, so such basins of attraction can be 
seen on discretizations. 

With what has just been said, we can already note that properties of discretizations of 
generic dissipative and conservative homeomorphisms are quite different: if before (almost) 
all possible dynamical properties are observed on the discretizations of a homeomorphism, 
the dissipative case reveals a much more regular behaviour: "the dynamics of discretizations 
of a homeomorphism / converges to the dynamics of /". For example, the measures /i-j^ 
accumulate on a single measure, namely /i;^, and not on all the invariant measures under 
/• 

10. The shredding lemma and his discrete counterpart 

In their article [AA12], F. Abdenur and M. Andersson try to identify the generic ergodic 
properties of continuous maps and homeomorphisms of compact manifolds. More precisely, 
they study the behaviour of Birkhoff limits for a generic homeomorphism / G Homeo(X) 
and almost every point x for Lebesgue measure. To do this they define some interesting 
behaviours of homeomorphisms related to Birkhoff limits, including one they call weird (see 
def 59). This definition is supported by their proof, based on the shredding lemma, that a 
generic homeomorphism is weird. We outline an improvement of this lemma, whose main 
consequence is that a generic homeomorphism has many open attractive sets, all of small 
measure and decomposable into a small number of small diameter open sets: 

Lemma 60 (Shredding lemma, F. Abdenur, M. Andersson, [AA12]). For all homeomor- 
phism f G Homeo(X)^ for all 5 > 0, there exists a family of regular pairwise disjoints open 
sets^^ /7i, . . . such that for all e' > 0, it exists g G Homeo(X) such that d{f,g) < S and: 

(t) giU]) C Uj, 
(it) \{Uj) < e, 

(in) \{\jW^j)>l-e, 
(xv) X{g{Uj))<sX{Uj), 

(v) there exists open sets Wj^i, . . . , Wj^e^ such that: 

(a) diam{Wj^i) < e' for all i e {1, • • . 

(b) g(W~i) C Wj, i+i, for every i € {1, . . . , £j - 1} and g(W~^) c Wj,i, 
(c) 

_ 

m>0 i=l 

The lemma as stated by F. Abdenur and M. Andersson in their paper can be obtained 
making e = e' . To get this improvement, it suffices to consider the sets Wj^i given by the 
weak version of the lemma and to crush them with a contraction mapping, as the authors 
do in the first part of the proof of the lemma. 

Remark 61. It can be easily seen that the lemma remains true on a neighbourhood of the 
homeomorphism g. 



*An open set is said regular if it is equal to the interior of its closure. 
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Remark 62. We can further assume that: 

(1) The open sets Wj^i have disjoints attractive sets, i.e. for ah j G N and ah i ^ i' ^ 
we have^^ 



(2) Each set Wj^i contains a periodic attractive point. Indeed, for ah e > 0, it exists 
X G Wj^i and m > such that d{x^f^{x)) < e. Then we perturb / in order to 
make x m-periodic and attractive, and this remains true in a neighborhood of /. 

(3) The set g{Uj) is "independent" from the choice of e': it exists compact sets Vj such 
that gijj)) C Vj and Vj C Uj, XiVj) < e\(Uj). 

This lemma tehs us a lot about the dynamics of a generic homeomorphism, which be- 
comes quite clear: there are many sinks whose basins of attraction are small and attract 
almost all the manifold X. Moreover there is convergence of the attractive sets of the 
shredding lemma to the closure of the sinks of /: 

Notation 63. For an homeomorphism / G Homeo(X), we note Aq the set of sinks of /, 
i.e. the sets of periodic points whose orbits attract nonempty open sets. 

Corollary 64. Let f G Homeo(X) verifying the conclusions of the shredding lemma for 
all £ = e' > {) and Wj^i^^ the corresponding open sets. Such homeomorphisms form a Gs 
dense subset o/Homeo(X). Then the sets 



converge for Hausdorff distance when e tends to to a closed set. This is a null set which 
coincides generically with the closure of the set Aq. 

Proof of corollary 64- Let / verifying the hypothesis of the corollary. We want to show 
that the sets A^ tend to Aq when e goes to 0. This is equivalent to show that for all (5 > 0, 
it exists £o > such that for all e < Sq, Aq C B{As^6) and A^ C B{Aq,S) (where B{A^S) 
denotes the set of points of X whose distance to A is smaller than S). Subsequently we will 
note Uj^s and Wij^s the open sets given by the shredding lemma for the parameters e = . 

We start by taking x G X whose orbit is periodic and attracts a nonempty open set, 
say O. Then there exists eo > such that for all e G]0^eo[^ there exists j G N such that 
the intersection between O and Uj^s is nonempty. Let y be an element of this intersection. 
Then there exists a subsequence of {f^{y))me'N which tends to x; moreover f^{y) G |Ji 
eventually. We deduce that x e [jj ■ VKj,i,£o all Sq small enough. Now, let (5 > 0. Since 
Aq is compact, it is covered by a finite number of balls of radius 6 centered in some points 
Xi whose orbits attract nonempty open sets. Taking the minimum of all the Eq associated 
to the Xi^ the inclusion Aq C B{A^^ 5) takes place for all s < £q. 

In the other way, let S > 0^ e < 5 and focus on the set By point 2. of the 

remark 62, we can suppose that there exists x G Wj^i^s whose orbit is periodic and attracts 
a nonempty open set. Thus x e Aq and since the diameter of is smaller than (5, 

Wj^i^s C B{Ao^S). The corollary is proved. □ 



These two sets are decreasing intersections of compact sets, so they are compact and nonempty. 
'And better, if we are given a countable family (Am)meN of good measures, generically Am(^o) = 0. 
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Now we establish a discrete counterpart of the shredding lemma. This result shows that 
the general dynamics of a generic dissipative homeomorphism is quite simple; thus the 
dynamics of the discretizations these homeomorphisms is also quite simple. More precisely, 
since having basins of attraction is stable by perturbation, we have a similar statement for 
discretizations of a generic homeomorphism. The following develops these arguments. 

To each point xn G E^^ we associate a closed set P^^{{xn})-, made of the points in X 
whose one of the projections on En is xn- The closed sets P^^{{xn}) form a basis of the 
topology of X when N runs through N and xn runs through En- Let / G Homeo(X), 

G N and : N ^ be a function that tends to +00 in +00. Let (5, e > and Uj be 
the sets obtained by the shredding lemma for /, 5 and e. 

For all j < £ we note Uj^ the union over xn e En of the closed sets P^^{{xn}) whose 
intersection with Uj are nontrivial. Then Uj^ tends to Uj for the metric d{A^ B) = X{AAB), 
and for the Hausdorff metric. By point (3) of remark 62, these convergences are independent 
from the choice of e' . Thus, for all k big enough, properties (i) to (iii) of the shredding 
lemma remain true for the discretizations gN (for arbitrary g satisfying the properties of the 
lemma). Taking e' small enough and modifying a little g if necessary, it exists sets Wj^i and 

g G Homeo(X) such that d{f,g) < 6, Wj^i C Pn^{{xn}) and Card(U^. ^ 1^^- ^ H^at) < i9{N). 
The others estimations over the sizes of the sets involved in the lemma are obtained similarly. 
Finally we have: 

Lemma 65 (Discrete shredding lemma). For all f G Homeo(X)^ for all e^S > and all 

function 1? : N ^ that tends to -\- 00 in -\-oo, it exists A/q G N such that for all N > Nq, 
it exists a family of subsets U^ ^ . . . , U^ of En and g G Homeo(X) such that d{f^g) < S 
and: 

(i) gN{Uf) c Uf, 
(it) Card(C/f ) < eqN, 

(iii) Card (U'=i f/f ) > (1 - s)qN, 

(iv) Card(5jv(C//)) <eCard([//), 

(v) for all j, it exists Wj'-^, . . . , Wj'f _ e E^- such that 

(a) Card(U, , < t?(iV)7 

(b) gN{wfi) = wj^i+i, for every i e {I, . . . ,tj - 1} and gN{wf^.) = wf^-^, 
(c) 



2>0 i=l 



(vi) for all j and all i, it exists sets Uj and Wj^i satisfying properties {i) to {v) of the 
shredding lemma such that Uj C P^"^(/7j^) and Wj^i C P^^(w^^^). 

(vii) if N > No, then for all j and all i, we have 



Y,dH{Uj,P^\Uf))<e and ^ ^i//(l^,,z, P^' (<,))< ^ 

3 hj 

( where dn is the Hausdorff metric ), and 

A(f/iAP^i(C/f )) < e and ^ \{W,,,AP^\wf^i)) < e. 

3 

Remark 66. Properties (i) to (v) are discrete counterparts of properties (i) to (v) of the 
continuous shredding lemma, but the two last ones reflect the convergence of the dynamics 
of discretizations to which of the original homeomorphism. 



This lemma implies that we can theoretically deduce the behaviour of a generic homeo- 
morphism from the dynamics of its discretizations. The next section details this remark. 
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11. Dynamics of discretizations of a generic homeomorphism 

To begin with, we deduce from the shredding lemma that the dynamics of discretiza- 
tions /at tends to that of the homeomorphism /. More precisely almost all orbits of the 
homeomorphism are (5-shadowed by the orbits of the corresponding discretizations. 

Definition 67. Let / and g be two maps from a metric space X into itself, x^y G X and 
^ > 0. We say that the orbit of x by / S -shadows the orbit of y by g if for all m G N, 
d(/™(x),5™(y))<,5. 

Corollary 68. For a generic homeomorphism f G Homeo(X)^ for all e > and all S > 0, 

it exists an open set A such that X{A) > 1 — e and A/q G N such that for all N > Nq and 
all X e A, the orbit of xn by f^ S-shadows that of x by f . 

Therefore, for a generic homeomorphism f , it exists a full measure dense open set O 
such that for all x G O, all 5 > and all N large enough, the orbit of xn by /at S-shadows 
that of X by f . 

Proof of corollary 68. This easily follows from the discrete shredding lemma, and especially 
from the fact that the points w^^^ tend to the sets Wj^i for Hausdorff metric. □ 

Therefore, we can "see" the dynamics of / on which of /at, with arbitrarily high precision 
provided that N is large enough. Among other things, this allows us to observe the basins 
of attraction of sinks of / on discretizations. Better yet, to each family of sinks {Wj^i)i 
of the basin Uj of the homeomorphism correspond a unique family of points {w^^^)i that 
are permuted cyclically by /at and which attracts a neighbourhood of Uj. Thus sinks are 
shadows by cyclic orbits of /at and we can detect the "period" of the sink (i.e. the integer 
kj) on discretizations. This behaviour is the opposite of what happens in the conservative 
case, where discretized orbits and true orbits are very different for most points. 

Again, in order to show that the dynamics of discretizations converge to that of the 
initial homeomorphism, we establish the convergence of attractive sets of /at to that of /. 
Recall that Aq is the closure of the sinks of / (see notation 63). 

Proposition 69. For a generic homeomorphism f G Homeo(X)^ the sets ^{fN) tend 
weakly to Aq in the following sense: for all e > 0, it existe Nq G'N such that for all N > Nq, 
it exists a subset En of En, stabilized by fN such that, noting 1^(/Ar) the corresponding 
maximal invariant set, we have card(^^) > 1 ~ ^ dni^o^ ^(/at)) < £■ 

Proof of proposition 69. Let £ > 0. For all G N, set En the union of the sets /7j^ of 
lemma 65 for the parameter e. This lemma ensures that En '^^ stable by /at and fills a 
proportion greater than 1 — e oi En- We also note ^^(/Ar) the associated maximal invariant 
set: 

^(/a^) = U ^/iv(^)- 

Property (viii) of lemma 65 ensures that 

]hR dH{A,,n{fN)) < e. 

To conclude, it suffices to observe that A^ Aq for Hausdorff distance. □ 

Finally, we conclude this section with a last consequence of lemma 65 which reflects 
that the ratio between the cardinal of the image of discretizations and which of the grid is 
smaller and smaller: 

Corollary 70. For a generic homeomorphism f G Homeo(X)^ 

^.^ Card(M^^Q. 

AT^+oo Card(£;Ar) ' 
particularly for all S > 0, it exists Nq e 'N such that if N > Nq, then the cardinal of each 
orbit of fN is smaller than 5Card(£^Ar). 
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This corollary can be seen as a discrete analogue of the fact that a generic homeomor- 
phism is totally singular, i.e. that there exists a Borel set of full measure whose image under 
/ is zero measure. Again, it reflects the regularity of the behaviour of the discretizations 
of a dissipative homeomorphism: generically, we can describe the behaviour of all (fine 
enough) discretizations. This is very different from the conservative case, where sometimes 
f^iE^) = En and sometimes Card(/Ar(£^Ar)) < 'd{N) where 'd : R+ is a given map 

that tends to +00 in +00. Moreover, we have the following estimation on the maximal 
invariant set (easily deduced from the discrete shredding lemma): 

Corollary 71. Let : R+ R^ be a function that tends to +00 in +00. Then for a 
generic homeomorphism f G Homeo(X)^ 

lim Cd^Td{n{fN)) <i^{N). 

It remains to study the behaviour of measures fi^^'^ (see definition 6). Again, the results 
are very different from the conservative case: for any open set /7, the measures ^ tend 
to a single measure, say /i^. 

Theorem 72. For a generic homeomorphism f G Homeo(X) and an open subset U of 
X, the measure is well defined^^ and is supported by the closure of the set of sinks Aq. 
Moreover the measures /j.'j^ jj tend weakly to . 

Sketch of proof of theorem 72. The proof of this theorem is based on the shredding lemma: 
the set of homeomorphisms which satisfy the conclusions of the lemma is a Gs dense, so 
it suffices to prove that such homeomorphisms / satisfy the conclusion of the proposition. 
Let U be an open subset of X and : X ^ R be a continuous function. 

We want to show that on one hand the integral (p d/i^ is well defined, i.e. that the 
Birkhoff limits for the function cp 

^ m—l 

are well defined for a.e. x G U; and on the other hand we have the convergence 

For the first step, the idea of the proof is that most of the points (for A) go in a Wj^i. 
Since the iterates of the sets Wj^i have small diameter, by uniform continuity, the function 
(p is almost constant on the sets f^{Wj^i). Thus the measure fil is well defined and almost 
constant on the set of points which goes in Wj^i. And by the same construction, since 
the dynamics of /at converge to that of /, and in particular that the sets Uj^ and {i^j^^} 
converge to the sets Uj and the measures ^ tend to the measures fi(j. □ 

Part 3. Numerical simulations 

12. Conservative homeomorphisms 

Now we present the results of numerical simulations of conservative homeomorphisms we 
conducted. We sought whether it is possible to observe behaviours such as those obtained 
in section 6 or that described by theorem 55 on actual simulations: it is not clear that the 
orders of discretization described by these results can be reached in practice, or if simple 
examples of homeomorphisms behave the same way as generic homeomorphisms. 

Recall that we simulate homeomorphisms f{x^y) = P o Q o P(x, ?/), where P and Q are 
both homeomorphisms of the torus that modify only one coordinate: 

P(x, y) = (x, y + p{x, y)) and Q{x, y) = {x ^ q{x, y),y), 



In other words a generic homeomorphism is weird, see also [AA12]. 
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and that we discretize according to the uniform grids on the torus: 

EN = [{^,...,'^)eT-\^j,0<ij<N-l 

To obtain conservative homeomorphisms, we choose p{x^y) = p{x) and q{x^y) = q{y). We 
perform simulations of two conservative homeomorphisms which are smah perturbations of 
identity or of the standard Anosov automorphism A : {x , y) {x -\- y , x -\- 2y): 

• To begin with we study fi = P o Q o P, with 

p{x) = cos(27r X 227x) + sin(27r x 253x), 

Q{y) = ^ cos(2^ X 241^) + ^ sin(2^ x 217^). 

This conservative homeomorphism is a smah perturbation of identity. Experi- 
ence shows that even dynamical systems with fairly simple definitions behave chaot- 
ically (see e.g. [GT83]). We can expect that a homeomorphism such as /i has a 
complex dynamical behaviour and even behaves essentially as a generic homeomor- 
phism, at least for small orders of discretization. Note that we choose coefficients 
that have (virtually) no common divisors to avoid arithmetical phenomena such as 
periodicity or resonance. 

• We also simulate /2 the composition of /i with the standard Anosov automorphism 
A, say /2 = PoQoPoA. This makes it a small C^-perturbation of A. Thus /2 
is semi-conjugated to A but not conjugated: to each periodic orbit of A corre- 
sponds many periodic orbits of /2. As for /i, we define /2 in the hope that the 
behaviour of its discretizations is fairly close to that of discretizations of a generic 
homeomorphism . 

Note that the homeomorphisms /i and /2 have at least one fixed point (for /i, simply 
make simultaneously p{x) and q{y) equal to 0; for /2, note that is a persistent fixed 
point of A). Thus, the theoretical results indicate that for a generic homeomorphism / 
which has a fixed point, infinitely many discretizations has a unique fixed point; moreover 
a subsequence of (/i;j^)Ar tends to an invariant measure under / supported by a fixed point. 
We can test if this is true on simulations. 

From a practical point of view, we restricted ourselves to grids of sizes smaller than 
2i5 ^ 2^^: the initial data become quickly very large, and the algorithm creates temporary 
variables that are of size of the order of five times the size of the initial data. For example, 
for a grid 2^^ x 2^^, the algorithm takes between 25 and 30 Go of RAM on the machine. 

12.1. Combinatorial behaviour. As a first step, we are interested in some quantities 
related to the combinatorial behaviour of discretizations of homeomorphisms. These quan- 
tities are: 

• the cardinal of the maximal invariant set Q{fN), 

• the number of periodic orbits of the map /at, 

• the maximal size of a periodic orbit of /n. 

Recall that according to corollaries 22, 27 and 39, for a generic homeomorphism, the ratio 
Card(Q(/jv)) j^^g^ fiuctuatc between and 1 depending on the number of orbits must 
fiuctuate between 1 and (e.g.) ^/qN and the maximal size of a periodic orbit must fiuctuate 
between 1 and qn- 

We calculated these quantities for discretizations of orders 128/c, for k from 1 to 100 and 
have represented them graphically. For information, if = 128 x 100, then qn — 1.6.10^. 

We begin with the cardinal of the maximal invariant set ^{/n) (figure 6). Contrary to 
what the theoretical results provide for a generic homeomorphism, for all simulations the 
ratio ^^^^g^^*^^^^ tends to when N increases. More specifically, the cardinal of ^{/n) 
evolves much more regularly for fi than for /2: for /i the value of this cardinal seems 
to be the sum of a smooth increasing function and a random noise, but for /2 this value 
seems to be the product of a smooth increasing function with a random noise. We have no 
explanation to the parabolic shape of the curve for /i: it refiects the fact that the cardinal 
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Figure 6. Size of the maximal invariant set Q{{fi)N) (left), number of periodic orbits of 
{fi)N (middle) and length of the largest periodic orbit of {fi)N (right) depending on 
for fi (top) and /2 (bottom), on the grids En with = 128/c, /c = 15, . . . , 100 



of l^((/i)Ar) evolves in the same way as -^/qn (whereas for a random map of a finite set 
with q elements into itself it evolves the same way as y^). Finally, it is interesting to note 
that for /2 , the size of the maximal invariant set is distributed more or less around the size 
of the maximal invariant set of a random map of a set with qn elements into itself, which 
depends (asymptotically) linearly of N and is worth about 16, 000 for N = 128 x 100 (see 
[BolOl]). 

According to the results of the previous sections, for a generic homeomorphism /, the 
number of periodic orbits of /at must fluctuate between 1 and (e.g.) y^^Ar- It is clearly 
not the case for the simulations. In addition, its behaviour is not the same for /i and 
for /2 (figure 6): for /i the number of orbits reaches quickly a value around 2.10^ to 
increase slightly thereafter, while for /2 it oscillates between 1 and 18, regardless the order 
of discretization. These graphics can be compared with those representing the size of the 
maximal invariant set Q{{fi)N)' the number of periodic orbits and the size of the maximal 
invariant set are of the same order of magnitude for /i (up to a factor 5), which means that 
the average period of a periodic orbit is small (which is not surprising, since /i is a small 
perturbation of identity). In contrast, they differ by a factor roughly equal to 10^ for /2, 
which means that this time the average period of a periodic orbit is very large. This can be 
explained partly by the fact that the standard Anosov automorphism tends to mix what 
happens in the neighborhood of identity. Fact remains that these simulations (such as the 
size of the maximal invariant set Q{{fi)N)) suggest that the behaviour in the neighborhood 
of identity and of the standard Anosov automorphism are quite different, at least for small 
orders of discretization. 

Regarding the maximum size of a periodic orbit of /n (figure 6), again its behaviour 
does not correspond to that of a generic homeomorphism: it should oscillate between 1 and 
qn^ which is not the case. However, it varies widely depending on especially when N 
is large. The qualitative behaviours are similar for both simulations, but there are some 
quantitative differences: the maximum of the maximal length of a periodic orbit is greater 
for /2 than for /i. 

12.2. Behaviour of invariant measures. We also simulated the measure for the two 
examples of conservative homeomorphisms /i and /2 as defined in page 31. Our purpose 
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is to test whether phenomena as described by theorem 55 can be observed in practice 
or not. We obviously can not expect to see the sequence (/i;j^)ArGN accumulating on all 
the invariant probability measures of /, since these measures generally form an infinite- 
dimensional convex set, but we can still test if it seems to converge or not. 

We present images of sizes 128 x 128 pixels representing in logarithmic scale the density 
of the measures /ij^: each pixel is coloured according to the measure carried by the set 
of points of £^ AT it covers. The blue corresponds to a pixel with very small measure and 
the red to a pixel with very high measure. Scales on the right of each image corresponds 
to the measure of one pixel in logarithmic scale to base 10: if the green codes —3, then 
a green pixel will have measure 10~^ for For information, when Lebesgue measure is 
represented, all the pixels have a value about —4.2. 
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Figure 7. Maximal value of the measure jj^n on Ejy depending on N for /i (left) and /2 
(right), on the grids En with N = 128^, = 15, . . . , 100 

Theoretically, the maximal value of the measure /j^n (figure 7) should oscillate between 
I/qn and 1, but this is not the case for these examples. Again the behaviour is quite 
different for /i and for for /i it seems that the values of jj^n are well distributed 
between null function and a linear function of N while for /2 we observe peaks: there 
are a few values of N for which jj^n is much higher than elsewhere. For /i the maximal 
value of /iAT seems globally grow with but with a quite irregular behaviour, apparently 
confirming theorem 55. Finally, we note that the maximum value of /in is much lower for 
/2 than for /i. 




The results of simulations of invariant measures of discretizations of /i (which is a 
conservative perturbation of identity) are quite positive: they agree with theoretical results 
about discretizations of generic conservative homeomorphisms, in particular with theorem 
55. When /i is discretized, we observe that the measure is first fairly well distributed. 
When the step of the discretization improves we can observe places where the measure ac- 
cumulates; moreover these places changes a lot when the orders of discretizations varies (see 
figure 8). This fairly agrees with what happens in the generic case, where the measure 
depends very much on the order of discretization rather than on the homeomorphism 
itself. There is also an other phenomenon: when the size of the grid is large enough (around 
10^^ X 10^^), it appears areas uniformly charged by the measure their sizes seems to 
be inversely proportional to the common mass of the pixels of the area. 




Figure 9. Simulations of invariant measures /i}^^ on grids of size N x with N = 
20005, . . . , 20013 (from left to right and top to bottom) 

For the discretizations of /2, which is a conservative perturbation of the standard 
Anosov automorphism, the simulations on grids of size 2^ x 2^ might suggest that the 
measures /i-j^ tend to Lebesgue measure. In fact, making a large number of simulations, 
we realize that there is also strong variations of the behaviour of measures (figure 9): the 
measure is often well distributed in the torus and sometimes quite singular with respect 
to Lebesgue measure (as in figure 7). For example when we discretize on the grid of size 
20010 X 20010 (middle right of figure 9), one observe an orbit with length 369 which mass 
84% of the total measure. In fact the behaviour of discretizations looks the same that in 
the neighbourhood of identity, modulo the fact that the standard Anosov automorphism 
tends to spread the attractive periodic orbits of discretizations on the entire torus: for many 
values of N composing by A spreads the behaviour of the measure , but sometimes (in 
fact seldom) a fixed point of {fi)N which attracts a large part of En is located around 
one of the few periodic points of small period for A. This then creates a periodic orbit for 
{f2)N with a big measure for 
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13. DiSSIPATIVE HOMEOMORPHISMS 

Let us finish by presenting ttie results of some numerical simulations of dissipative tiomeo- 
morptiisms. Again, our aim is to compare the theoretical results with the reality of numer- 
ical simulations: for simple homeomorphisms and reasonable orders of discretization does 
one have convergence of the dynamics of discretizations to that of the homeomorphism, as 
suggested by the above theorems? 

Recall that we simulate homeomorphisms such that 

/(x, y) = Qo P(x, y) or f{x, y) = P o Q o P{x, y), 

where P and Q are two homeomorphisms of the torus that modify only one coordinate: 

P(x, y) = (x, y + p{x, y)) and Q{x, y) = {x ^ q{x, y),y). 

Unlike the conservative case, p and q depends on both x and y. Again, we tested two 
homeomorphisms : 

• To begin with we studied fs = Q o P^ with 

p{x, y) = ^ cos(2^ X 227y) + ^ sin(2^ x 233x), 

q{x, y) = ^ cos(2^ X 241^) + ^ sin(2^ x 217^) + ^ cos(2^ x 271x) 

This dissipative homeomorphism is a small perturbation of identity, whose de- 
rivative has many oscillations whose amplitudes are close to 1. That creates many 
fixed points which are sinks, sources or saddles. 

• It has also seemed useful to us to simulate a homeomorphism close to the identity 
in topology, but with a small number of sinks. Indeed, as explained heuristically 
by J.-M. Gambaudo and C. Tresser in [GT83], a homeomorphism such that fs can 
have a large number of sinks whose basins of attraction are small. It turns out that 
the dissipative behaviour of fs can not be detected in orders discretization that can 
be achieved in practice. We therefore defined another homeomorphism close to the 
identity in topology, but with much less sinks, say = P o Q o P^ with 

p{x, y) = (^0 cos(27r X ^)) + ^ (50 cos(27r x 5x)) , 

q{x, ^) = ^ th (50 cos(2^ x ^)) + ^ th (50 cos(2^ x 7^)) + ^ th (50 cos(2^ x 3x)) . 

13.1. Combinatorial behaviour. We simulated some quantities related to the combina- 
torial behaviour of discretizations of homeomorphisms, namely: 

• the cardinal of the maximal invariant set ^^(/tv), 

• the number of periodic orbits of /at, 

• the maximal size of a periodic orbit of /at. 

We calculated these quantities for discretizations of orders 128/c for k from 1 to 100 and 
represented it graphically. For information, if TV = 128 x 100, then qn — 1.6.10^. 

Theoretically, the ratio between the cardinal of ^^(/Ar) and qn must tend to 0, this is 
what can be seen in simulations. This is not really surprising: we already observed that 
for discretizations of conservative homeomorphisms. In this context it is interesting to 
compare the behaviour of ^^(/Ar) in the conservative and the dissipative case. The result 
is a bit disappointing: for /a, the graphics in the case of conservative and dissipative 
homeomorphisms are as alike as two peas in a pod, while in theory they should be very 
different. Still, for /4, there are very few sinks, and around each one there are just a few 
attractive points of the discretization (each sink is "sharp") so that the cardinal of l^((/4)Ar) 
is more or less constant. For this case, the behaviour of discretizations seems clearly typical 
of the dissipative case. 

The behaviour of the number of periodic orbits of /at is not provided by the theoretical 
study. However, we can hope that it reflects the fact that the dynamics converges to that 
of the initial homeomorphism: among others, we can test if it is of the same order as 
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Figure 10. Size of the maximal invariant set Q{{fi)N) (left), number of periodic orbits 
of {fi)N (middle) and length of the largest periodic orbit of {fi)N (right) depending on 
for fs (top) and (bottom), on the grids En with N = 128/c, /c = 15, . . . , 100 



the number of sinks of the homeomorphism. In practice, its behaviour is similar to that 
of the cardinal of the maximal invariant set (^(/at): its evolution is very regular for fs 
but oscillates (from = 15) between 30 and 90 for f^. As in the conservative case, the 
comparison with figure 10 indicates that the average length of a periodic cycle is about 3 
for fs and f^. 



Since the dynamics of discretizations is assumed to converge to that of the initial homeo- 
morphism, one could expect that the length of the longest periodic orbit of discretizations 
{fi)N is almost always a multiple of that of an attractive periodic orbit of fi. The graphic 
of this length for fs looks like the conservative case, so one can say that the dissipative 
behaviour of this homeomorphism is not detected for reasonable orders of discretizations. 
Despite this, the values of the lengths of the longest periodic orbit are much smaller (up 
to a factor 10) than in the conservative case. It is not obvious if it is more a dissipative 
effect than a coincidence. For F4, the length of the longest orbit is smaller than 25, which 
is normal because we are supposed to detect the sinks of the homeomorphism. The changes 
of the length of the longest orbit are probably due to errors of discretization if the neigh- 
borhood of the sinks of /4. Furthermore, we note that some values of the size of the longest 
orbit are more achieved: for example for 30 different integers N the longest periodic orbit 
has length 4, for 19 different integers N this length is 8 and for 10 different integers it 
is 10. This is certainly due to the shadowing property of periodic orbits of f^ by those of 
the discretizations. 



13.2. Behaviour of invariant measures. As in the conservative case we calculated the 
invariant measures /i-j^ of dissipative homeomorphisms fs and f^ as defined on page 36. 
Our aim is to test whether theorem 72 apply in practice or if there are technical constraints 
such that this behaviour can not be observed on these examples. For a presentation of the 
representations of the measures, see page 33. 



38 



PIERRE-ANTOINE GUIHENEUF 




10 20 30 40 50 GO 70 80 90 100 10 20 30 40 50 GO 70 80 90 100 



Figure 11. Maximal value of /i at on En depending on N for fs (left) and (right), on 
the grids En with N = 12Sk, = 15, . . . , 100 

There is no clear theoretical results concerning the behaviour of the maximal value of 
the measure /ij^ (figure 11), but we can at least observe that a convergence of this amount 
would be a good indication that the considered homeomorphism is dissipative. As in 
the conservative case (in the neighbourhood of identity), the behaviour of /i^^ seems very 
chaotic. While is always smaller than 0.07, the maximum value of fi-j^ is very high and 
fluctuates between 0.02 and 0.3, with many values around 0.06 and 0.12 and some values 
around 0.25. The obvious arithmetic relationships between these values can be explained 
by the fact that some sinks of discretizations merge for some values of N. Thus, there 
are orbits which attract many points; it seems to reflect the typical behaviour of generic 
homeomorphisms . 




The behaviour of invariant measures for /a, which is a smah dissipative perturbation 
of identity, is relatively similar to that of invariant measures for fi i.e. the corresponding 
conservative case: when the order discretization is large enough, there is a strong varia- 
tion of the measure /i^^. Moreover this measure has a significant absolutely continuous 
component with respect to Lebesgue measure. Nevertheless, there are differences with the 
conservative case: the maximum value of /ij^ is very large, much more than for /i, and 
there are orbits which attract many points. On the other hand, one can observe the follow- 
ing phenomenon (we hope that it is anecdotic): these orbits are located on vertical stripes. 
This is probably related to specific arithmetic phenomena due to the specific form of the 
homeomorphism /a. 

It is normal not to be able to observe all the sinks of fs on these simulations: as noted by 
J.-M. Gambaudo and C. Tresser in [GT83], the size of these sinks can be very small com- 
pared to the numbers involved in the definition of fs. So even in orders discretization such 
as 2^^, the dissipative nature of the homeomorphism can not be detected on discretizations. 
For this reason, it looks as if the discretizations of fs are very similar to those of /i. 




Figure 13. Simulations af /n^j^ on the grids En, with N = 2^ and /c = 6, . . . , 14 (from left 
to right and top to bottom) 

Recall that what happens for /s is pretty close to what happens for /i. For their part, 
simulations of invariant measures for on grids of size 2^ x 2^ highlight that we expect from 
a generic dissipative homeomorphism: the measures /i^^ tend quickly to a single measure 
(that is also observed on a series of simulations), which is carried by the sinks of /4. The 
fact that /4 has few sinks allows the discretizations of reasonable orders (typically 2^^) to 
find the actual sinks of the initial homeomorphism, contrary to what we had observed for 

/3. 
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